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ABSTRACT 


An  analysis  of  the  integral  wavelength  direct  normal 
incidence  solar  radiation  and  its  spectral  components  passing 
through  the  RG2  and  RG8  filters  has  been  made  for  measurements 
taken  at  Edmonton,  Alberta,  Canada  during  the  period  May  1972 
to  September  1975  under  clear  sky  conditions.  The  basic  terms 
and  equations  useful  in  solar  energy  calculations  are  presented. 
The  annual  variations  of  the  apparent  solar  constant  and  the 
apparent  extinction  coefficient  have  been  examined.  The 
diurnal  and  annual  variations  of  the  direct  radiation  and  of 
its  spectral  components  have  been  studied.  The  average  spec¬ 
tral  composition  of  the  direct  radiation  is  found  to  be  56% 
in  the  infrared  region,  37.5%  in  the  shortwave  region  and  6.5% 
in  the  red  region  for  solar  altitudes  between  60  and  20  degrees. 
Expressing  the  infrared  and  shortwave  fractions  as  functions 
of  the  optical  air  mass  resulted  in  linear  relationships.  For 
the  same  altitude,  the  atmospheric  transmission  factors  for 
the  direct  radiation  and  its  spectral  components  are  lower  in 
summer  than  in  winter.  The  average  ratios  of  the  extinction 
coefficients  for  the  infrared  and  the  shortwave  regions  to  that 
for  the  integral  wavelength  are  0.79  and  1.15  during  the 
summer,  respectively.  The  corresponding  values  in  the  winter 
are  0.66  and  1.4.  The  atmospheric  turbidity  and  the  water 
vapor  absorption  have  been  determined  and  they  show  marked 
annual  variations  with  maximum  values  during  the  summer. 
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CHAPTER  1 


PRINCIPLES  OF  SOLAR  ENERGY 
CALCULATIONS 


1 
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1.1  INTRODUCTION 

The  sun  is  one  of  the  main  sources  of  natural  energy, 
and  since  this  energy  source  is  undepleted  while  most  of  the 
other  sources  of  energy  are  depleting,  the  utilization  of 
solar  energy  in  as  many  everyday  living  applications  as 
possible  is  desirable.  Engineering  and  scientific  problems 
involved  with  utilizing  solar  energy  on  the  earth  must  be 
studied,  and  in  many  countries  research  on  fundamental  and 
applied  topics  on  harnessing  of  solar  energy  in  thermal, 
electrical  and  chemical  processes  is  expanding. 

In  solar  radiation  problems  a  quantitive  knowledge  of 
incident  solar  energy  is  needed  for  all  types  of  solar  weather 
including  clear,  partly  cloudy,  and  cloudy  days.  One  of  the 
prime  requirements  is  the  prediction  of  the  annual  variation 
of  solar  energy  on  an  oriented  surface  at  any  location  on 
the  earth's  surface.  A  surface  receives  the  maximum  solar 
energy  under  clear  sky  conditions.  Thus,  an  investigation 
was  undertaken  to  measure  and  correlate  the  annual  variation 
in  the  direct  beam  solar  energy  under  clear  sky  conditions 
that  will  be  received  on  a  surface  oriented  perpendicular  to 
the  sun's  rays  at  Edmonton, Alberta .  In  addition  the  distri¬ 
bution  of  energy  within  specific  wavelength  bands  of  the 
solar  spectrum  was  also  studied. 

In  sections  1.2  to  1.5  of  this  chapter,  a  brief  review 
is  given  about  the  nature  of  the  sun  (the  energy  source)  and 
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the  earth  (the  receiver) .  Also  the  basic  relationships  and 
essential  terms  required  in  the  calculations  are  defined. 

In  section  1.6  of  this  chapter  the  procedure  for  collecting 
data,  the  basic  calculations  and  specific  information  about 
the  selected  location  are  presented.  In  chapter  2  the  trans¬ 
mission  of  the  direct  solar  beam  through  the  earth's  atmos¬ 
phere  and  the  energy  carried  within  the  whole  spectrum  will 
be  studied.  Chapter  3  will  deal  with  the  filter  measurements 
of  solar  radiation  and  the  energy  distribution  in  specific 
wavelength  regions  reaching  the  ground.  Finally,  in  chapter 
4  the  turbidity  of  the  atmosphere  and  the  energy  absorbed  by 
the  atmospheric  water  vapor  will  be  determined  from  the  solar 
radiation  measurements.  Conclusions  are  summarized  in 
chapter  5. 

1.2  THE  SUN  AND  THE  EARTH 

1.2.1  THE  SUN:  One  of  the  stars  of  the  universe.  The 
diameter  of  the  sharp  circular  boundary  of  the  sun  is  about 
865,400  miles,  and  its  mass  is  about  332,000  times  that  of 
the  earth.  It  is  entirely  gaseous.  Astrophysicists  generally 
divide  the  structure  of  the  sun  into  three  regions: 
i)  the  interior  which  is  the  main  mass  of  the  sun  is  believed 
to  have  a  central  temperature  of  many  millions  of  degrees 
and  a  central  gas  pressure  of  perhaps  a  billion  atmospheres. 
It  is  here  that  the  sun's  energy  output  is  generated.  Astro¬ 
physicists  have  established  that  the  sun's  main  constituent 
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is  hydrogen  and  it  is  believed  that  the  sun's  energy  output 
occurs  from  conversion  of  hydrogen  into  helium  in  the 
presence  of  carbon  and  nitrogen. 

ii)  The  photosphere  is  a  thin  layer  of  gas  which  forms  the 
bright  boundary  of  the  sun.  Both  pressure  and  density  are 
very  low  in  this  layer  and  the  temperature  is  about  10/000°F. 
The  greatest  part  of  the  sun's  thermal  radiation  is  emitted 
from  the  photosphere. 

iii)  The  solar  atmosphere  is  composed  of  the  chromosphere  and 
corona.  Both  regions  of  the  solar  atmosphere  offer  little 
resistance  to  the  sun's  radiation  output.  The  chromosphere 

is  a  relatively  thin  layer  of  gas  at  extremely  low  temperature 
and  low  density.  The  corona  which  is  the  extremity  of  the 
sun,  consists  of  rarified  gases  which  extend  beyond  the 
chromosphere  for  a  million  miles  or  more. 

1.2.2  THE  EARTH:  The  planet  earth  is  nearly  spherical 
with  a  diameter  of  about  7,900  miles.  It  makes  one  rotation 
about  its  axis  every  24  hours  and  completes  an  orbit  about 
the  sun  in  a  period  of  365.25  days  approximately.  The  earth's 
orbit  about  the  sun  is  approximately  circular,  with  the  sun 
located  slightly  off  the  center  of  the  circle.  The  earth's 
mean  distance  from  the  sun  is  about  92,900,000  miles  and 
about  January  1,  the  earth  is  closest  to  the  sun  while  about 
July  1,  it  is  most  remote,  the  variation  being  3.3  percent. 

The  earth's  axis  of  rotation  is  tilted  23.5  degrees  with 
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respect  to  its  orbit  plane  and  is  of  profound  significance, 
for  together  with  the  earth's  daily  rotation  and  yearly 
revolution  it  accounts  for  the  distribution  of  solar  radia¬ 
tion  over  the  earth's  surface,  the  changing  length  of  hours 
of  daylight  and  darkness,  and  the  changing  of  the  seasons. 


1.3  BASIC  EARTH-SUN  ANGLES 

The  zenith  angle,  altitude  angle,  azimuth  angle  and 
incidence  angle  are  earth-sun  angles  used  in  solar  calcula¬ 
tions.  They  may  be  calculated  for  any  location  P  on  the 
earth's  surface  if  the  latitude  L  and  hour  angle  h  for  that 
location,  and  the  sun's  declination  6  are  known.  It  is  con¬ 
venient  to  illustrate  these  angles  using  the  astronomical 
coordinate  systems  [27]. 

The  systems  of  astronomical  spherical  coordinates  are 
all  analogous  to  the  system  of  latitudes  and  longitudes  on 
the  earth.  The  most  important  difference  being  that  the 
earth  is  viewed  from  the  outside  and  the  celestial  sphere 
which  is  a  fictitous  sphere  of  indeterminate  radius,  with 
its  center  at  the  observer,  and  the  sun  is  imagined  to  be  pro¬ 
jected  on  the  inside  surface  of  it  is  viewed  from  the  inside. 
Every  such  system  contains  two  diametrically  opposite  poles, 
a  north  pole  and  a  south  pole,  the  north  pole  being  the  one 
that  is  visible  from  the  north  pole  of  the  earth.  There  is 
a  primary  circle,  which  is  a  great  circle  of  the  sphere,  so 
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placed  that  each  pole  is  90  degrees  from  it.  The  latitudinal 
coordinate  of  an  object  is  its  angular  distance  measured 
northward  or  southward  from  the  primary  circle.  The  longi¬ 
tudinal  coordinate  of  an  object  is  the  angular  distance 
measured  along  the  primary  circle  from  a  specified  origin 
to  the  point  of  intersection  with  a  great  circle  passing 
through  the  object  and  the  two  poles.  Three  basic  systems, 
the  terrestrial  system,  and  the  two  astronomical  systems 
most  used  for  locating  the  sun  (the  equatorial  system  and 
the  horizon  system),  are  compared  in  table  1.1.  The  equa¬ 
torial  system.  Fig.  1.1,  is  the  geometrical  extension  of 
the  terrestrial  system.  The  horizon  system  shown  in  Fig.  1.2 
is  fixed  relative  to  the  observer. 

Mow  it  is  easy  to  define  the  three  basic  angles,  the 
latitude  L,  the  local  hour  angle  h,  and  the  sun's  declination 
<$ .  Figure  1.3  illustrates  these  angles. 

1.3.1  LATITUDE, L:  The  angular  distance  of  the  point 
P  north  or  south  of  the  earth's  equator. 

1.3.2  HOUR  ANGLE,  h:  The  angle  through  which  the  earth 
must  turn  to  bring  the  meridian  of  P  directly  in  line  with 
the  sun's  rays.  At  solar  noon  the  hour  angle  is  zero.  The 
hour  angle  is  measured  on  either  side  of  solar  noon  and  one 
hour  of  time  equals  360/24  or  15  degrees  of  earth's  rotation. 
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TABLE  1 . 1  Comparison  of  the  terrestrial 
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The  terrestrial  meridian  is  that  half  of  a  great  circle  on  earth  that  passes  through  points 
having  the  same  longitude  and  terminates  at  the  north  and  south  poles. 


TABLE  1.1 
(continued) 


(2)  The  celestial  meridian  is  a  great  circle  on  the  celes¬ 
tial  sphere  passing  through  the  two  celestial  poles 
and  the  observer's  zenith. 


(3)  The  zenith  is  the  point  directly  overhead. 


(4)  The  altitudinal  coordinate  is  latitude  if  it  is  measured 
from  the  horizon  or  zenith  if  it  is  measured  downward 
from  the  zenith. 


declination  of 


Meridian  of  observer 


north 

celestial  pole 


segment  of  hour 
through  the  sun 


N 


local  hour  angle  of  the  sun 


celestial  equator 
(prime  circle) 


horizon  circle 


Figure  1 . 1  Equatorial  astronomical  coordinate  system. 
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segment  of  vertical  circle 


Figure  1 .2  Horizon  system  of  astronomical  coordinates. 


Figure  1.3  Illustration  of  the  three  basic  earth-sun  angles; 
Latitude  L,  Hour  angle  h  and  Decl ination  %  . 
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1.3.3  DECLINATION,  6:  The  sun's  declination  is  the 
angular  distance  of  the  sun  north  or  south  of  the  celestial 
equator.  The  tilt  of  the  earth's  axis  is  the  direct  cause 
of  the  sun  declination.  The  sun's  declination  would  be  23.5 
degrees  south  of  the  equator  at  the  winter  solstice,  23.5 
degrees  north  of  the  equator  at  the  summer  solstice,  and  at 
the  equinoxes,  the  declination  is  zero.  Because  the  period 
of  the  earth's  complete  revolution  about  the  sun  does  not 
coincide  exactly  with  a  calendar  year,  the  declination  varies 
slightly  for  the  same  day  from  year  to  year.  A  declination 
formula  is  given  by  [14]: 

6  =  23.45  sin  [  (360/365)  (284  +  n) ] 

where  n  is  the  day  of  the  year  measured  from  January  1.  The 
sun's  declination  may  be  obtained  from  Solar  Ephemeris  Publica¬ 
tions  [22 ]  . 

1.3.4  DERIVED  SOLAR  ANGLES:  With  the  sun's  declination 
6,  the  hour  angle  h,  and  the  latitude  L  for  a  specific  locality 
known,  the  sun's  zenith  z,  altitude  a,  azimth  angle  <j>  and 
incidence  angle  0  can  be  calculated  using  the  following  equa¬ 
tions  : 


cos  z  =  cos  L  •  cos  h  •  cos  6  +  sin  L  •  sin  6  (1.1) 


Since  the  altitude  angle  a  =  (tt/2)  -  the  zenith  angle  z 


therefore,  sin  a  =  cos  L  •  cos  h  •  cos  6  +  sin  L  •  sin  5  (1.2) 
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and  the  azimuth  angle  cp ,  measured  from  south  is  given  by- 

cos  cf>  =  sec  a  (cos  L  •  sin  6  -  cos  5  •  sin  L  •  cos  h) 

(1.3) 

In  applying  equations  1.1  to  1.3  north  latitudes  are  con¬ 
sidered  positive  and  south  latitudes  negative.  The  declina¬ 
tion  will  be  positive  for  the  summer  period  between  the 
vernal  equinox  and  autumnal  equinox  (March  22  to  September  22 
approximately)  and  negative  at  other  times.  The  hour  angle 
in  degrees  is  measured  on  either  side  of  solar  noon.  Thus 
h  is  limited  to  values  between  zero  and  tt.  If  h  <  tt/2, 
cos  h  is  positive  and  if  h  >  tt/2,  cos  h  is  negative. 

For  a  particular  surface  orientation  the  sun's  incidence 
angle  9  is  given  by: 

cos  0  =  cos  i  •  sin  a  +  sin  i  •  cos  a  •  cos(c}>  -  s) 

where:  i  is  the  inclination  angle  of  the  surface  under 

consideration,  and  is  expressed  in  degrees  from 
the  horizontal,  and 

s  is  the  aspect  of  the  surface,  and  is  expressed  in 
degrees  from  the  south  with  a  westerly  component 
being  positive. 

For  example  a  vertical  east  facing  surface  would  have  an 
inclination  angle  of  90  degrees  and  an  aspect  of  -90  degrees. 
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The  solar  angles  z,  a,  <j> ,  the  incidence  angle  9  and 
the  surface  angles  i  and  s  can  be  defined  through  the  geomet¬ 
rical  relationships  between  the  sun  and  the  horizontal  plane 
and  the  inclined  surface  at  the  point  of  observation  on  the 
ground.  Figure  1.4  illustrates  these  angles.  The  three 
reference  axes  are  Z (vertical),  W(west)  and  S (south)  and  P 
is  the  point  of  observation,  then  by  definition 

Line  AP  is  the  sun's  ray 

Line  BP  is  the  projection  of  the  sun's  ray  on  the 
horizontal  plane 

Line  CP  is  the  normal  to  the  inclined  surface 

Line  DP  is  the  projection  of  the  normal  to  the  inclined 
surface  on  the  horizontal  plane. 

Sun  altitude  angle,  a:  the  angle  in  the  vertical  plane 
between  the  sun's  ray  and  the  projection  of  the  sun's  ray  on 
the  horizontal  plane  (angle  APB) , 

Sun  zenith  angle  z:  the  angle  between  the  sun's  rays 
and  a  line  perpendicular  to  the  horizontal  plane  (angle  APZ) , 

Sun  azimuth  angle,  cf> :  the  angle  in  the  horizontal  plane 
measured  from  south  to  the  horizontal  projection  of  the  sun's 
rays,  in  the  northern  hemisphere,  (angle  BPS) , 

Incident  angle,  0:  the  angle  between  the  sun's  rays  and 
the  normal  to  the  surface  (angle  APC) , 

Surface  inclination  angle,  i,  is  measured  from  the  hori¬ 


zontal  plane  . 


13 


Figure  1 .4  Illustration  of  solar  angles. 


. 
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Surface  azimuth  angle,  s  (also  known  as  the  aspect  of 
the  surface) :  It  is  the  angle  in  the  horizontal  plane 
measured  from  south  to  the  horizontal  projection  of  the  normal 
line  to  the  surface  (angle  DPS) . 


1.4  SOLAR  TIME 

Solar  radiation  calculations  must  be  made  in  terms  of 
solar  time,  usually  with  respect  to  solar  noon  which  occurs 
with  the  sun  at  the  zenith  for  the  station  meridian.  In 
order  to  calculate  mean  solar  noon  from  true  solar  noon  some 
terms  should  be  clearly  defined. 

1.4.1  TRUE  SOLAR  TIME,  TST :  The  time  that  would  be 
shown  by  a  sun  dial,  or  it  is  the  hour  angle  of  the  sun  and 
it  is  measured  by  the  diurnal  motion  of  the  sun.  An  apparent 
solar  day  at  any  place  is  the  interval  of  time  between  two 
successive  transits  of  the  sun  over  the  meridian  of  that 
place,  with  the  transit  from  midnight  to  midnight  as  being 
more  adaptable  to  the  civil  day.  True  solar  time  is  the 
known  angle  of  the  sun  plus  12  hours.  Whereas  a  civil  day 
is  precisely  24  hours,  a  true  solar  day  is  slightly  different 
due  to  irregularities  of  the  earth's  rotation,  obliquity  of 
the  earth's  orbit  and  other  factors. 

1.4.2  MEAN  SOLAR  TIME,  MST:  It  is  measured  by  the 
diurnal  motion  of  a  fictitious  body  called  the  mean  sun,  which 
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is  supposed  to  move  uniformly  in  the  celestial  equator, 
completing  the  circuit  in  one  tropical  year.  A  mean  solar 
day  at  any  place  is  the  interval  of  time  between  two  succes¬ 
sive  transits  of  the  mean  sun  over  the  meridian  of  that 
place.  Since  mean  solar  time  is  uniform  and  regular  in  its 
passage,  clocks  and  watches  are  regulated  to  it.  Mean  solar 
time  is  called  civil  time. 

1.4.3  THE  EQUATION  OF  TIME,  ET:  The  equation  of  time 
is  the  difference  in  hour  angle  between  true  sun  and  mean 
sun,  or  it  is  the  quantity  added  to  mean  solar  time  to  obtain 
true  solar  time.  The  true  sun  is  sometimes  ahead  and  some¬ 
times  behind  the  mean  sun  by  an  amount  which  varies  from 
zero  to  about  16  minutes. 

TST  =  MST  +  ET  (1.4) 

1.4.4  GREENWICH  CIVIL  TIME,  GCT:  Time  reckoned  from 
midnight  at  the  greenwich  meridian  (zero  longitude)  is  known 
as  Greenwich  civil  time  or  universal  time.  Such  time  is 
expressed  on  an  hour  scale  from  zero  to  24.  Thus,  midnight 
is  0  and  noon  is  12. 

1.4.5  LOCAL  CIVIL  TIME,  LCT:  The  mean  solar  time  rec¬ 
koned  from  the  precise  longitude  of  the  observer.  At  the 
same  instant,  LCT  is  more  advanced  by  four  minutes  for  each 
degree  difference  in  longitude  east  of  Greenwich  meridian  and 
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is  less  advanced  by  four  minutes  for  each  degree  difference 
in  longitude  west  of  Greenwich. 

1.4.6  ZONE  AND  STANDARD  TIMES ,  ZST:  To  avoid  the 
inconvenience  of  the  continuous  change  of  mean  solar  time 
with  longitude,  zone  time  or  civil  time  is  used.  The  earth 
is  divided  into  24  zones  each  approximately  15  degrees  wide, 
and  centered  on  standard  longitudes  0,  15,  30  and  so  on  until 
180  degrees  east  and  west  of  Greenwich.  Within  each  of  these 
zones  the  clock  time  is  the  same  for  all  locations  of  the 
zone  and  is  kept  at  the  mean  solar  time  of  the  standard 
meridian. 

Now  the  true  solar  time  TST  at  any  instant  at  any  loca¬ 
tion  on  the  earth's  surface  of  longitude  G  can  be  computed 
from  equation  1.4,  substituting  for  MST  as  follows: 

TST  =  MST  +  ET 
=  LCT  +  ET 

=  ZST  +  (GS  -  G)  *  4  +  ET  (1.5) 

where  GS  is  the  longitude  of  the  standard  meridian  of  the  zone 
including  that  location.  The  values  of  ET  can  be  obtained 
from  Nautical  Almanac  or  Solar  Ephemeris  publications. 
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1-5  THE  SOLAR  CONSTANT ,  THE  SOLAR  SPECTRAL  IRRADIANCE  AND 

THE  AIR  MASS: 

1.5.1  THE  SOLAR  CONSTANT,  1^:  The  total  solar  energy 
received  per  unit  time  per  unit  area  exposed  normally  to  the 
sun's  rays  at  the  mean  sun-earth  distance  in  the  absence  of 
the  terrestrial  atmosphere. 

The  value  of  the  solar  constant  which  had  been  proposed 
by  F.S.  Johnson  in  1954  [21]  is  442  Rtu  ft~2hr_1  (2  cal  cm"2min_1)  . 
Recently  a  new  value  of  429  Btu  ft-2  hr"1  (1.94  cal  cm"2  hr"1) 
has  been  proposed  for  the  solar  constant  [33] .  This  is  3% 
lower  than  the  Johnson  value. 

1.5.2  SOLAR  SPECTRAL  IRRADIANCE:  The  distribution  of 
the  solar  constant  as  a  function  of  wavelength.  The  spec¬ 
trum  of  the  sun  extends  from  x-rays  of  wavelength  1  Angstrom 
or  less  to  radio  emission  of  wavelength  100  meters  and  more. 
However,  the  wavelength  range  0.275  to  4.95ym  contains  99% 
of  the  energy  and  99.9%  is  within  the  range  0.220  to  11.8ym. 

The  proposed  standard  for  zero  air  mass  solar  spectral  irradi- 
ance  which  is  presented  in  reference  [33]  was  used  in  this 
research. 

1.5.3  THE  AIR  MASS,  m:  The  air  mass  is  the  ratio  of  the 
length  of  the  path  of  the  sun's  rays  through  the  atmosphere 

to  the  length  of  the  path  when  the  sun  is  at  the  zenith.  In 
solar  radiation  calculations,  unit  depth  of  the  atmosphere  is 


' 
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taken  as  the  depth  when  the  sun  is  at  the  zenith  for  sea 
level . 

For  solar  altitude  angles  greater  than  or  equal  to  20 
degrees,  the  air  mass  is  equal  to  the  cosecant  of  the  altitude 
angle  a: 


Air  mass  =  1/sin  a,  a  -  20  degrees  (1.6) 

Below  20  degrees  a  more  accurate  method  should  be  used  and 
the  formula  proposed  by  Robinson  for  calculating  the  air  mass 
of  a  homogeneous  spherical  atmosphere  is  often  used  [27]. 

It  is  given  by: 

2  ^ 

Air  mass  =  [  (r/H  •  sin  a)  +2  (r/H  )+l]  -  r/H  •  sin  a  (1.7) 

Where  r  is  the  radius  of  the  earth,  H  =  P  /p  g 

o  o  ^ 

Where  P  is  the  pressure  of  the  earth's  surface  (sea  level), 
o 

and  p  is  the  density  at  the  earth's  surface.  On  substituting 
o 

Pq  =  760  mm  Hg  and  Pq  corresponding  to  the  temperature  of  273°K, 
it  is  found  that  H  =7.991  km. 

The  air  mass  calculated  from  equation  1.6  or  1.7  is  the 
relative  optical  air  mass  m^.  The  absolute  optical  air  mass  m 
can  then  be  computed  using  the  relation: 

m  =  mh  •  P/pQ  (1.8) 

where  p  is  the  surface  pressure  at  the  place  and  time  of  obser¬ 


vation  . 
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1.6  DATA  AND  COMPUTATIONS 

The  data  in  this  study  are  the  values  of  the  normal  in¬ 
cidence  insolation  and  the  insolation  passing  the  RG2  and 
RG8  filters  measured  by  Eppley  normal  incidence  pyrheliometers 
under  clear  sky  conditions.  The  Pyrheliometers  are  mounted 
on  the  roof  of  the  Mechanical  Engineering  building  of  the 
University  of  Alberta  campus  at  Edmonton,  Alberta  and  the 
output  is  recorded  by  a  strip  chart  recorder.  A  continuous 
record  is  available  during  the  months  of  April  to  September, 
but  only  a  partial  record  is  available  during  the  months  of 
October  to  March.  Clear  sky  conditions  during  both  forenoon 
and  afternoon  are  difficult  to  obtain  in  the  region  as  there 
is  a  tendency  for  a  clear  sky  during  the  forenoon,  but  light 
fluffy  clouds  often  appear  just  before  solar  noon  and  these 
may  remain  for  the  rest  of  the  day  [28] .  About  200  days  were 
chosen  from  a  4-year  period,  (1972  to  1975) ,  with  about  70% 
of  them  during  the  months  of  May  to  September.  The  days  selec¬ 
ted  showed  a  high  value  of  normal  incident  insolation  at  mid¬ 
day,  and  a  symmetrical  shape  for  forenoon  and  afternoon. 
Whenever  possible,  full  day  data  was  used  in  the  analysis  but, 
due  to  infrequency  of  occurrence ,  these  days  were  supplemented 
by  half-day  data,  and  part-day  data  that  were  deemed  to  be 
clear  sky  conditions. 

The  sun  crosses  the  local  meridian  at  local  apparent 
noon  each  day,  i.e.  the  local  hour  angle  is  zero  and  the  true 
solar  time  is  12.  Edmonton  (113°  31'W  and  53°  34'N)  is 
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included  in  the  zone  with  a  standard  meridian  of  105°W  and  the 
local  standard  time,  LST,  corresponding  to  the  true  Solar  noon 
on  any  given  day  was  computed  from  equation  1.8  as  follows: 

12  =  LST  +  [4*  (105°  -  113°  31'  )  +  ET]/60 

LST  =  12  +  (34. 08  -  ET)/60  (1.9) 

This  means  that  on  any  given  day  when  the  sun  crosses  the 
local  meridian  at  Edmonton,  the  local  clock  indicates  12  hour 
and  34.08  -  ET  minutes.  In  other  words,  the  true  solar  time 
lags  the  local  standard  time  by  34.08  -  ET  minutes.  The 
values  of  the  equation  of  time  ET  substituted  in  equation  1.9 
were  obtained  from  Solar  Ephemeris  [22]  and  the  chart  records 
were  adjusted  to  true  solar  time.  The  readings  were  taken 
from  the  charts  every  15  minutes  whenever  it  was  possible, 
starting  from  solar  noon  as  hour  angle  zero  till  hour  angle 
6  for  both  forenoon  and  afternoon  periods  during  the  months 
of  April  to  September  and  to  hour  angle  4  during  the  rest  of 
the  year.  These  values  were  entered  on  computer  cards  and 
the  quantities  of  the  solar  energy  were  computed  for  normal 
incidence  insolation  and  insolation  passing  the  RG2  and  RG8 
filters.  It  was  found  that  the  data  showed  only  a  small 
change  during  any  15-minute  interval,  so,  only  the  values  at 
the  exact  hour  angle  from  0  to  6  were  used  in  most  of  the 
analysis.  The  remaining  data  was  utilized  to  augment  the 
material  for  certain  calculations  especially  when  the  best 
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fit  was  obtained  by  the  computer. 

For  each  day  the  altitude  of  the  sun  was  calculated 
every  15  minutes  from  equation  1.2  using  the  declination  values 
published  in  Solar  Ephemeris  [22].  The  corresponding  relative 
optical  air  mass  (termed  hereinafter  as  relative  air  mass) 
was  calculated  using  equations  1.6  and  1.7,  then  corrected  to 
absolute  optical  air  mass  (termed  hereinafter  as  optical  air 
mass)  according  to  equation  1.8.  For  this  correction  the 
average  daily  pressure  values  for  forenoon  and  afternoon  were 
obtained  from  Upper  Air  Bulletins  [11] .  In  these  Bulletins 
pressure  values  are  published  for  Stony  Plain,  Alberta  (776 
meters  above  mean  sea  level)  and  a  correction  for  elevation 
level  was  applied  to  get  the  pressure  at  Edmonton  (668  meters 
above  mean  sea  level) . 

For  some  calculations,  it  was  necessary  to  reduce  the 

measured  values  of  the  solar  insolation  to  the  mean  distance 

Rq  between  the  sun  and  the  earth  by  multiplication  by  the 

2 

factor  (R^/Rq)  ,  obtained  from  the  tables  published  in  Nautical 

Almanac  [25],  where  R^  is  the  earth-sun  distance  at  time  of 


observation. 


. 


CHAPTER  2 


CLEAR  SKY  DIRECT  NORMAL  INCIDENCE 
SOLAR  RADIATION 
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2 . 1  SOLAR  EXTRATERRESTRIAL  RADIATION 

The  extraterrestrial  solar  radiation  at  normal  incidence 
at  any  instant,  1^,  will  differ  from  the  solar  constant,  1^, 
because  of  the  variation  in  the  earth-to-sun  distance 
throughout  the  year.  This  results  in  a  maximum  variation  of 
irradiation  of  ±3.5  percent.  To  obtain  1^  the  solar  constant 
should  be  multiplied  by  the  factor  (RQ/R^)  .  On  January  1 

the  value  of  is  7  percent  higher  than  on  July  1. 

The  solar  radiation  is  depleted  as  it  penetrates  the 
earth’s  atmosphere.  The  intensity  of  the  direct  beam  received 
at  the  earth's  surface  depends  upon  the  path  length  traversed 
by  the  sun  rays.  The  length  of  the  path  through  the  atmos¬ 
phere  depends  on  the  solar  declination,  the  time  of  the  day, 
and  the  location  on  the  earth's  surface.  The  extent  to  which 
any  particular  air  mass  depletes  the  sun  rays  depends  on  the 
components  of  the  earth's  atmosphere. 


2 . 2  ATMOSPHERIC  EFFECTS  ON  EXTRATERRESTRIAL  ENERGY 

Scientific  studies  of  the  effects  of  the  earth's  atmos¬ 
phere  upon  solar  radiation  has  shown  that,  in  general,  the 
atmospheric  elements  absorb  and  scatter  part  of  the  incident 
solar  energy.  Scattering  may  be  accomplished  by  dry  gas 
molecules  (nitrogen,  oxygen,  and  other  molecular  components 
of  the  atmosphere),  by  dust,  aerosols  and  water  droplets. 
Small  amounts  of  many  substances  absorb  solar  energy,  but  the 
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principal  absorbing  mediums  are  ozone,  water  vapor,  carbon 
dioxide,  and  cloud  particles. 

The  remaining  portion  of  the  original  direct  radiation 
may  reach  the  earth's  surface  unchanged  in  wavelength  and  is 
known  as  the  direct  solar  radiation  or  normal  incidence  insola¬ 
tion.  A  portion  of  the  scattered  energy  will  return  to  space, 
but  some  of  it  also  reaches  the  earth's  surface  as  diffused 
radiation.  It  comes  from  the  entire  sky  vault  and  therefore 
is  usually  termed  diffuse  sky  radiation.  On  clear  days 
(cloudless  sky) ,  diffuse  sky  radiation  is  small  compared  to 
direct  solar  radiation  but  during  extremely  cloudy  days,  only 
diffuse  radiation  may  reach  the  ground.  In  general  a  surface 
on  the  earth  receives  solar  energy  associated  with  radiation 
of  two  forms;  direct  radiation  and  diffuse  radiation.  This 
research  is  concerned  with  the  incidence  of  direct  solar 
radiation  on  clear  days  (cloudless  sky) . 


2. 3  EXTINCTION  OF  SOLAR  RADIATION  AND  ATMOSPHERIC 
TRANSMISSION  FACTOR 

The  intensity  I ,  of  parallel  monochromatic  energy  trans- 

A 

mitted  through  the  earth's  cloudless  atmosphere  is  given  by  [27] 


I 


(2.1) 


where  A  is  the  monochromatic  extinction  coefficient, 

A 


m  is  the  optical  air  mass. 
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^OA  '*'S  intensity  of  parallel  monochromatic  beam  of 
wavelength  A  outside  the  atmosphere,  air  mass  zero, 

I-v  is  the  intensity  of  the  parallel  beam  at  the  ground. 


The  atmospheric  transmission  T-^  of  a  parallel  mono¬ 
chromatic  beam  of  wavelength  A  through  the  atmosphere  gives 
the  fraction  of  energy  incident  outside  the  earth's  atmos¬ 
phere  (air  mass  zero)  which  penetrates  the  atmosphere  to  the 
earth's  surface  along  optical  air  mass  m.  From  equation  2.1 
it  is  given  by: 


T 


A 


(2.2) 


(2.3) 


where  is  the  monochromatic  transmittance  with  the  sun  in 
the  zenith  i-e  for  unit  air  mass. 

Several  investigators  have  treated  the  monochromatic 
extinction  coefficient  in  different  ways,  but  as  a  com¬ 
promise,  it  includes  the  effect  of  both  scattering  and 
absorption.  They  all  have  considered  scattering  due  to  main 
scattering  agents  in  the  atmosphere,  dry  gas  molecules,  dust 
and  water  vapor.  Also  they  agreed  that  absorption  is  mainly 
by  ozone,  and  water  vapor  in  absence  of  clouds.  On  these 
bases  the  extinction  coefficient  A,  can  be  written  to 
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represent  the  sura  of  four  parameters  as  follows: 


AX  =  aRU)  +  aDU)  +  aW(A)  +  a0(A)  (2.4) 

where  a^(A)is  the  extinction  coefficient  for  clean  dry  air 
and  it  accounts  for  scattering  by  air  molecules 
according  to  Rayleigh's  theory, 

a^(A)is  the  extinction  coefficient  of  the  atmospheric 
dust.  It  represents  extinction  due  to  turbidity 
of  the  atmosphere,  taking  into  account  scattering 
and  absorption  due  to  dust  as  well  as  scattering 
by  water  vapor  in  the  atmosphere, 

a^(A)is  the  coefficient  of  absorption  by  the  water 
vapor  contained  in  the  atmosphere,  and 

aQ(A)is  the  extinction  coefficient  due  to  ozone. 

There  is  rarely  more  than  0.4  centimeters  of  ozone  in  the 
vertical  column  of  the  atmosphere  and  it  is  concentrated 
mainly  at  elevations  between  15  and  35  kilometers  above  the 
ground.  Yet  the  absorption  coefficient  of  ozone  is  so  large 
that  the  spectrum  of  solar  energy  at  the  ground  is  cut  off 
below  0.29  micron.  Therefore,  from  the  view  point  of  utili¬ 
zation  of  solar  energy  near  the  ground  our  attention  is  found 
on  wavelengths  longer  than  0.29  micron.  Then  equation  2.4 
can  be  written  as: 


aR(X)  +  aD(X)  +  a^A) 


(2.5) 
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The  calculation  of  the  extinction  coefficients  in  equation 
2.4  is  illustrated  in  reference  [12].  The  extinction  due  to 
atmospheric  turbidity,  dust,  and  absorption  by  water  vapor 
in  the  atmosphere  will  be  discussed  in  chapter  4. 

If  equation  2.1  is  integrated  over  all  the  solar  spect- 

2 

rum  and  the  factor  (R^/R^)  is  used,  the  basic  relation  for 
the  extinction  of  solar  radiation  through  the  earth's  atmos¬ 
phere  may  be  written  as: 


I 


<VRi> 


-A ( X ) •m 


VA) 


dX 


0 


(2.6) 


where  Rq  is  the  mean  earth-sun  distance, 

R.  is  the  earth-sun  distance  at  time  of  observation, 

l 

I  ( X )  is  the  solar  constant  as  a  function  of  the 
wavelength  X , 

I  is  the  intensity  of  the  direct  solar  radiation  at 
the  earth's  surface, 

m  is  the  optical  air  mass,  and 

A ( X )  is  the  extinction  coefficient  and  it  is  composed 
of  four  components  depending  on  the  wavelength  X 
as  it  is  given  in  equation  2.4. 


2.3.1  TRANSMISSION  COEFFICIENT  OR  ATMOSPHERIC 

TRANSMISSION  FACTOR  FOR  DIRECT  SOLAR  RADIATION: 

This  is  defined  as  the  fraction  of  the  energy  outside  the 
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earth's  atmosphere  that  is  received  at  the  earth's  surface 
as  normal  incidence  radiation.  Therefore  it  is  the  result 
of  the  integration  of  equation  2.6  at  certain  air  mass  divided 
by  the  extraterrestrial  radiation  at  air  mass  zero. 

2 . 4  APPARENT  SOLAR  INSOLATION  AT  AIR  MASS  ZERO  AND  APPARENT 

ATMOSPHERIC  EXTINCTION  COEFFICIENT 

Although  equation  2.1  is  theoretically  applicable  only 
for  monochromatic  radiation,  it  is  suggested  to  be  applicable 
for  the  total  direct  solar  radiation  in  the  method  proposed 
by  the  American  Society  of  Heating,  Refrigerating,  and  Air 
Conditioning  Engineers  [  1  ]  of  expressing  solar  energy  quan¬ 
tities.  This  method  represents  the  solar  insolation  at  the 
earth's  surface,  under  clear  sky  conditions,  as  a  function 


* 


of  the  apparent  solar  insolation  at  air  mass  zero,  I  ,  and 
the  atmospheric  extinction  coefficient,  A.  The  normal  direct 
incidence  insolation  is  represented  by  the  equation: 


I 


(2.7) 


where  I  is  irradiation  on  a  plane  normal  to  the  solar 


beam  at  the  ground. 


* 

I  is  the  apparent  solar  irradiation  at  air  mass  zero, 
or  the  apparent  solar  constant. 


A 


is  the  apparent  extinction  coefficient,  and 
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m  is  the  optical  air  mass. 

The  logarithm  of  equation  2.7  yields  a  straight  line 
relationship  between  I  and  m  and  when  plotted  on  semi  loga- 
rithmic  paper,  has  an  intercept  equal  to  I  and  the  slope 
equal  to  A. 

This  method  implies  that  the  atmospheric  extinction 

coefficient  does  not  vary  over  the  wavelength  and  also  is 

constant  during  the  period  of  measurements.  To  the  extent 

that  this  is  not  practically  true,  equation  2.1  when  applied 

to  experimental  data  of  the  total  direct  solar  spectrum, 

* 

results  in  an  ordinate  I  which  is  less  than  the  solar  con¬ 
stant  1^.  The  error  introduced  by  applying  this  method  is 
difficult  to  estimate.  However,  if  the  extinction  coefficient 
does  not  have  a  large  variation  over  the  wavelength  regions 
that  contain  the  major  portion  of  the  solar  energy,  equation 
2.7  will  give  results  that  can  be  used  with  sufficient  accu- 

•k 

racy  for  most  design  situations  [28].  The  values  of  I  and 
A  in  equation  2.7  vary  during  the  year  because  of  seasonal 
changes  in  the  dust  and  water  vapor  content  of  the  atmosphere, 
and  also  because  of  the  changing  earth-sun  distance. 

The  high  altitude  solar  irradiance  measurement  from  an 
air  craft  [32]  indicate  that  the  plot  of  the  logarithm  of 
normal  incidence  versus  air  mass  produces  a  curve  which  is 
concaved  up  rather  than  a  straight  line.  Measurements  at 
ground  level  tend  to  support  this  curvature.  Although  the 
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curvature  is  small,  representing  the  curve  by  a  best  fit 

straight  line  over  a  certain  air  mass  range  will  give  a 

* 

different  value  of  I  and  A  than  if  the  fit  is  taken  over  a 

different  air  mass  range.  Thus  in  using  the  straight  line 

* 

extrapolation  to  zero  air  mass  method  the  value  of  I  and  A 
is  influenced  by  the  air  mass  range  used  in  the  computation 
and  the  local  atmospheric  conditions.  Consideration  should 
be  given  to  the  air  mass  range  used  in  the  calculation  when 
comparing  different  stations'  results  or  even  different 
results  at  the  same  station. 


2 . 5  MEASUREMENT  OF  THE  APPARENT  SOLAR  CONSTANT  AND  THE 

APPARENT  EXTINCTION  COEFFICIENT  AT  EDMONTON,  ALBERTA, 

CANADA 

Taking  the  natural  logarithm  of  equation  2.7  results  in 
the  straight  line  relationship 

In  (I)  =  -A* *  m  +  In ( I*)  (2.8) 

which  has  a  negative  slope  of  magnitude  A,  and  an  ordinate 

•k 

value  of  In (I  )  at  air  mass  zero. 

The  annual  variation  at  Edmonton,  of  apparent  solar  in- 

* 

solation  at  air  mass  zero,  I  ,  and  the  apparent  extinction 
coefficient.  A,  were  calculated  using  equation  2.8.  About 
80  of  the  clearest  days  in  the  period  between  May  1972  and 
September  1975  were  selected  for  analysis.  The  criteria  for 
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selection  were:  (1)  high  I  value  at  solar  noon,  (2)  close 
symmetrical  shape  for  morning  and  afternoon  and  (3)  full 
day  data  supplemented  by  half  day  data.  The  days  that  were 
selected  represent  all  seasons,  but  there  were  more  days  in 
the  months  May  to  September  than  the  remaining  months.  This 
may  be  a  result  of  the  short-term  record  or  it  may  indicate 
that  the  frequency  of  clear  sky  is  higher  in  the  period  May 
to  September.  For  each  selected  day  the  measured  values  of 
the  normal  incidence  insolation,  IDN,  were  prepared  at  hour- 
angle  0,  1  and  at  successive  15-minute  intervals  thereafter 
till  hour-angle  6  for  the  months  of  April  to  September  and  to 
hour-angle  4  for  the  months  of  October  to  March.  The  corres¬ 
ponding  values  of  the  optical  air  mass  were  computed  as 
described  in  section  1.6.  About  16  to  22  data  points  were 
available  for  each  day  (for  each  forenoon  and  afternoon) 
during  the  months  of  April  to  September,  but  only  8  to  16 
data  points  for  each  day  during  the  other  months.  The  air 
mass  and  insolation  values  were  entered  on  cards  and  the  best 
straight  line  for  equation  2.8  was  fitted  by  computer,  using 
a  least  squares  technique.  A  separate  line  was  calculated 
for  each  forenoon  and  afternoon  and  the  calculations  were 
cut  off  at  air  mass  8. 

The  resulting  values  for  the  apparent  extinction  coeffi¬ 
cient  A,  and  the  apparent  solar  insolation  at  air  mass  zero 
* 

I  are  presented  in  figures  2.1  and  2.2  respectively.  The 
solid  curves  were  fitted  by  eye  through  the  data  points  to 


■< 
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FIGURE  2.1  ANNUAL  VARIATION  OF  THE  APPARENT  EXTINCTION 


COEFFICIENT  .A.  MEASURED  AT  EDMONTON* 


' 
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FIGURE1  2.2  ANNUAL  VARIATION  OF  THE  APPARENT  SOLAR  INSOLATION 


AT  AIR  HASS  ZERO 


HEASURED  AT  EDHONTON • 
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show  the  annual  variation  of  both  A  and  I  . 

A  similar  study  was  done  by  Sadler  [28]  at  the  same  loca¬ 
tion  using  data  for  the  period  1965-1968.  He  represented  the 

annual  variation  of  the  apparent  extinction  coefficient.  A, 

* 

and  the  apparent  solar  insolation  at  air  mass  zero,  I  ,  by 
the  following  two  periodic  curves,  respectively: 

A  =  0.134  +  0.056  cos  [360  (n  -  155) /365 ]  (2.9) 

I*=  351.4  +  21.0  cos  [360  (n  -  130)/365]  (2.10) 

where  n  =  day  of  year  measured  from  1  January.  These  periodic 

curves  are  shown  as  broken  lines  on  figures  2.1  and  2.2 

* 

Table  2.1  lists  the  average  values  for  A  and  I  obtained 
from  the  periodic  curves  and  from  the  present  study  for  each 
month  of  the  year.  In  general,  a  good  agreement  seems  to 
exist  between  the  two  studies.  However,  the  present  results 

•k 

show  about  5%  lower  values  for  I  for  the  months  April  to 
August  inclusive.  A  better  agreement  appears  for  the  rest  of 
the  year.  Comparing  values  for  A,  the  present  study  shows 
slightly  lower  values  for  January  through  September  and 
slightly  higher  values  for  October  through  December. 


I 


. 
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TABLE  2 . 1  Comparison  between  Sadler's  and  present  monthly 

* 

average  values  for  A  and  I  obtained  at  Edmonton. 


Month 

A 

I*,  Btu.ft  2. hr  1. 

Sadler 

Present 

Sadler 

Present 

January 

0.096 

0.085 

345 

342 

February 

0.123 

0.103 

356 

356 

March 

0.149 

0.125 

365 

363 

April 

0.174 

0.152 

371 

353 

May 

0.188 

0.164 

372 

350 

June 

0.188 

0.171 

367 

347 

July 

0.172 

0.152 

358 

340 

August 

0.146 

0.141 

347 

330 

September 

0.117 

0.110 

337 

335 

October 

0.092 

0.100 

331 

334 

November 

0.079 

0.087 

331 

334 

December 

0.080 

0.083 

335 

336 

2. 6  ANNUAL  VARIATION  OF  THE  DIRECT  SOLAR  RADIATION  AT 
NORMAL  INCIDENCE ,  IDN',  AT  VARIOUS  HOUR  ANGLES  AS 
MEASURED  AT  EDMONTON 


The  serial  hourly  values  of  the  clear  sky  direct  normal 
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incidence,  IDN,  solar  radiation  measured  at  Edmonton  were 

made  available  in  computer  form  as  described  in  section  1.6. 

—2 

Figures  2.3  to  2.9  show  these  values  in  Btu*ft  *hr  for  hour 
angles  zero  to  six  throughout  the  year.  The  curves  in  these 
figures  were  best  fitted  by  eye  to  show  the  approximate  annual 
variation  of  IDN.  Figure  2.10  depicts  all  seven  curves  on  one 
graph  for  visual  comparison.  The  annual  variation  at  all  hour 
angles  shown  is  partly  because  of  the  variation  of  the  earth- 
sun  distance  throughout  the  year,  and  because  of  seasonal 
change  in  the  dust  and  water  vapor  content  of  the  atmosphere. 
Dust  and  water  vapor  are  the  major  factors  attenuating  solar 
radiation  on  its  path  through  the  atmosphere. 

In  Figure  2.10,  at  hour  angle  zero  (solar  noon)  the 
intensity  of  the  direct  normal  incidence  radiation  shows  a 
slight  variation  from  March  through  September.  The  maximum 
occurs  in  March  and  April.  The  lower  values  of  IDN  from 
October  to  February  are  primarily  due  to  the  longer  path 
through  the  atmosphere.  The  curves  at  hour  angles  one  and 
two  are  similar  and  close  to  the  solar  noon  curve  because 
of  the  small  change  in  the  path  lengths  for  hour  angles  zero, 
one  and  two.  At  hour  angles  3,  4,  5  and  6  the  maximum  energy 
is  available  during  May  to  August.  The  available  energy 
decreases  rapidly  in  the  other  months  due  to  rapid  increases 
in  the  path  length  of  the  direct  radiation  through  the 
atmosphere.  The  difference  between  the  energy  measured  at 
two  consecutive  hours  increases  after  hour  angle  two  for  any 


-  * 


IDN  (BTU./SQ  FOOT/HOUR  ) 

140.00  180.00  220.00  260.00  300.00  340.00 
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FIGURE  2.3  ANNUAL  VARIATION  OF  THE  CLEAR  SKY 

DIRECT  NORMAL  INCIDENCE  SOLAR 
RADIATION  .IDN.  MEASURED  AT 


EDMONTON  AT  HOUR  ANGLE  ZERO  . 
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IDN  (BTU./SQ  FOOT/HOUR  1 

140.00  180.00  220.00  260.00  300.00  340.00 
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FIGURE  2-4  RNNUflL  VRRIRTI0N  OF  THE  CLERR  SKY 

DIRECT  NORMAL  INCIDENCE  SOLRR 
RADIATION  .IDN.  MEASURED  AT 


EDMONTON  AT  HOUR  ANCLE  1  . 
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IDN  (BTU./SQ  FOOT/HOUR  ) 

140.00  180.00  220.00  260.00  300.00  340.00 
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FIGURE  2.5  ANNURL  VARIATION  OF  THE  CLEAR  SKY 

DIRECT  NORMAL  INCIDENCE  SOLAR 
RADIATION  .IDN,  MEASURED  AT 


EDMONTON  AT  HOUR  ANGLE  2  . 


IDN  (BTU./SQ  FOOT/HOUR  ) 

140.00  180.00  220.00  260.00  300.00  340.00 
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FIGURE  2.6  ANNUAL  VARIATION  OF  THE  CLEAR  SKY 

DIRECT  NORMAL  INCIDENCE  SOLAR 
RADIATION  .IDN.  MEASURED  AT 


EDMONTON  AT  HOUR  ANGLE  3  . 


IDN  (BTU./SQ  FOOT/HOUR  ) 

140.00  180. 00  220.00  260.00  300.00  340.00 
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FIGURE  2.7  RNNURL  VRRIRTION  OF  THE  CLERR  SKY 

DIRECT  NORMRL  INCIDENCE  SOLRR 
RRDIRTION  .IDN.  NERSURED  RT 
EDMONTON  RT  HOUR  RNGLE  4-  . 


I DN  (BTU./SQ  FOOT/HOUR  ) 

140.00  180.00  220.00  260.00  300.00  340.00 
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FIGURE  2-S  ANNURL  VARIATION  OF  THE  CLEAR  SKY 

DIRECT  NORMAL  INCIDENCE  SOLAR 
RADIATION  *  I DN  *  MEASURED  AT 


EDMONTON  AT  HOUR  ANGLE  5  . 
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IDN  ( BTU  ./SQ  FOOT/HOUR  ) 

140.00  180.00  220.00  260.00  300.00  340.00 
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FIGURE  2-9  ANNUAL  VARIATION  OF  THE  CLEAR  SKY 

DIRECT  NORMAL  INCIDENCE  SOLAR 
RADIATION  .IDN.  MEASURED  AT 


EDMONTON  AT  HOUR  ANGLE  <o  . 


.  . 


IDN  (BTU./SQ  FOOT/HOUR  ) 

140.00  180.00  220.00  260.00  300.00  340.00 
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FIGURE  2-10  RNNURL  VARIATION  OF  THE  CLEAR  SKY 

DIRECT  NORMAL  INCIDENCE  SOLAR 
RADIATION  .IDN.  MEASURED  AT 
EDMONTON  AT  HOUR  ANGLES  O  to  G . 


particular  day.  This  is  mainly  due  to  the  rapid  change  of 
the  sun's  altitude  which  results  in  a  large  change  of  the 
air  mass  (path  length  through  the  atmosphere)  in  a  short 
period  of  time. 

Table  2.2  lists  the  monthly  average  of  direct  normal 

incidence,  IDN,  radiation  measured  at  Edmonton  at  various 

hour  angles.  These  values  occur  at  about  the  third  week  of 

-2  -1 

the  month.  Individual  days  may  differ  by ±  15  Btu*ft  *hr 
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TABLE  2 » 2  Monthly  average  normal  incidence  solar  radiation 

(IDN)  measured  at  Edmonton  at  different  hour 

-2  -1 

angles.  Unit:  Btu.ft.  hr. 


Month 

Hour 

Angle 

From  Solar  Noon 

0 

1 

2 

3 

4 

5 

6 

January 

255 

250 

227 

180 

February 

284 

279 

264 

233 

189 

March 

291 

286 

275 

262 

232 

178 

April 

294 

291 

284 

272 

250 

216 

152 

May 

291 

286 

285 

278 

262 

236 

194 

June 

290 

286 

285 

277 

265 

241 

210 

July 

288 

284 

282 

274 

258 

238 

204 

August 

286 

284 

279 

270 

250 

222 

168 

September 

286 

285 

280 

269 

243 

200 

153 

October 

281 

272 

262 

243 

201 

154 

November 

252 

251 

235 

193 

170 

December 

234 

228 

198 

156 

* 


CHAPTER  3 


DISTRIBUTION  OF  SOLAR  ENERGY  IN 
DIFFERENT  WAVELENGTH  REGIONS 


47 
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3. 1  THE  SOLAR  RADIATION  AT  THE  EARTH'S  SURFACE 

Although  the  spectrum  emitted  by  the  sun  extends  from 

o  °  -10 

x-rays  of  wavelength  of  a  fraction  of  an  Angstrom  (1A  =  10 

meter)  to  radio  emission  of  hundreds  of  meters  wavelength, 

only  the  spectral  region  from  0.3  micron  to  3.0  microns 
_  £ 

(ly  =  10  meter)  reaches  the  earth's  surface  as  direct  solar 
radiation.  Almost  all  the  energy  associated  with  radiation 
of  wavelength  less  than  0.24y  is  absorbed  by  oxygen  in  the 
upper  atmosphere  and  this  leads  to  the  formation  of  ozone. 

At  about  60  kilometers  above  the  ground  the  atmosphere 
becomes  so  rich  in  ozone  that  all  the  energy  in  the  wavelengths 
shorter  than  0.3y  is  absorbed  and  ground  measurements  are  not 
possible.  The  emitted  solar  energy  at  wavelengths  less  than 
0 . 3y  is  0.485  percent.  At  wavelengths  longer  than  3.0y  in 
the  infrared  range,  which  accounts  for  2.17  percent  of  the 
emitted  solar  energy,  most  of  the  energy  is  absorbed  by  atmos¬ 
pheric  water  vapor  and  carbon  dioxide. 

As  mentioned  in  chapter  two  the  intensity  of  the  direct 
normal  solar  radiation  incident  on  a  given  surface  on  the  earth 
depends  on  the  sun's  elevation  and  shows  a  diurnal  and  annual 
variation.  It  also  depends  on  its  attenuation  along  its  path 
through  the  earth's  atmosphere.  The  attenuation  of  the  solar 
beam  by  the  atmosphere  is  not  uniform  over  the  solar  spectrum. 
The  influences  which  cause  diminution  in  the  incident  energy 
passing  through  the  atmosphere  may  be  divided  into  two  major 


effects : 


' 

' 
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(i)  An  effect  which  is  generally  ascribed  to  the  scattering 
and  general  absorption  in  the  atmosphere.  It  affects  all 
wavelengths,  and  in  fact  the  shortwaves  more  than  the  long 
ones.  It  is  due  to  two  factors.  One  is  the  scattering  by  air 
molecules  (Rayleigh  scattering) .  The  other  factor  is  that 
part  of  the  incident  energy  scattered  and  absorbed  by  the  dust 
in  the  atmosphere. 

(ii)  The  second  effect  is  limited  to  certain  wavelengths  of 
the  spectrum  and  is  called  the  selective  absorption  by  the 
gases  of  the  atmosphere.  This  absorption  is  practically 
negligible  for  the  wavelength  range  0.3  to  0 . 6  y  and  the  largest 
absorption  occurs  in  the  infrared  region  (X  >  0.63y)  of 
the  solar  spectrum  and  it  is  mainly  due  to  water  vapor. 

Between  0.3  and  3.0  microns,  considerable  variation  in 
attenuation  of  the  solar  beam  occurs  for  the  different  wave¬ 
lengths.  Therefore,  it  is  desirable  to  divide  the  solar  spec¬ 
trum  into  several  wavelength  regions  and  study  the  distribu¬ 
tion  of  energy  in  each  region  in  addition  to  the  total  energy 
that  reaches  the  ground. 


3. 2  CHARACTERISTICS  OF  CUTOFF  GLASS  FILTERS  USED  IN  SOLAR 
RADIATION  MEASUREMENTS 


For  the  separation  of  the  solar  spectrum  into  reasonably 
well  defined  bands,  filters  are  more  practical  and  easier  to 
handle  than  a  spectrophotometer  [9].  Absorbing  glass  filters 


v 

■ 


* 
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can  be  readily  employed  with  normal  incidence  pyrheliometers 
by  simple  positioning  of  the  filter  in  front  of  the  detector, 
and  the  transmitted  radiation  can  be  expressed  directly  in 
absolute  units.  The  filters  recommended  for  use  by  the  Inter¬ 
national  Radiation  Commission  are  the  three  Schott  glass  filters 
OG1,  RG2  and  RG8 . 

Each  of  the  Schott  glass  filters  OG1,  RG2  and  RG8  has  a 
main  transmission  band  in  which  an  ideal  filter  shows  almost 
uniform  transparency  and  shows  no  transmission  outside  this 
region.  Practically,  each  filter  has  a  sharp  lower  cutoff 
separating  the  wavelengths  almost  totally  transmitted  from 
those  totaly  absorbed.  The  transmission  change  occurs  over 
a  20  to  50  my  range  and  the  center  of  the  lower  cutoff  is 
defined  by  the  wavelength  A  at  which  the  transmission  of  the 
filter  is  50  percent.  The  center  of  cutoff  for  filters  with 
the  same  main  transmission  characteristics  but  made  of  different 
melts  may  differ  in  the  position.  Also  it  is  slightly  tempera¬ 
ture  dependent. 

The  transmission  characteristics  of  each  filter  may  be 
expressed  as  a  filter  factor  and,  by  multiplying  the  measured 
radiation  behind  the  filter  by  this  factor,  the  solar  radiation 
in  front  of  the  filter  is  obtained.  If  the  filter  is  used 
with  a  pyrheliometer  having  a  weather  protective  window  in 
front  of  the  receiver,  the  filter  factor  has  to  be  reduced  by 
one  percent  to  compensate  for  the  effect  of  the  multiple 
reflections  between  the  filter  and  the  window. 


V 


. 


51 


To  compare  measurements  taken  at  different  locations  or 
taken  by  filters  made  of  different  melts,  a  correction  may 
be  made  to  account  for  the  shift  in  the  positions  of  the 
center  of  the  lower  cutoff  for  the  filter  from  the  cutoff 
position  of  a  standard  filter.  This  correction  is  made  to 
reduce  the  measurement  to  that  which  would  be  obtained  by  the 
standard  filter.  One  method  of  applying  this  correction  is 
to  test  the  filter  in  use  by  a  national  standarizing  labora¬ 
tory  to  obtain  a  correction  to  the  filter  factor.  This  method 
is  applicable  for  shifts  of  less  than  10  my  in  the  position  of 
the  lower  cutoff,  and  may  usually  be  applied  to  filters  of  the 
same  type  but  from  different  melts  [10].  Another  method,  which 
calculates  the  correction  as  a  function  of  the  shift  in  the 
center  of  cutoff,  the  air  mass,  and  the  turbidity  coefficient 

O 

of  the  atmosphere,  has  been  introduced  by  Angstrom  and 
Drummond  [9] .  They  suggested  that  if  the  shift  is  more  than 
10  my,  the  correction  is  evidently  so  large  that  it  ought  to 
be  taken  into  account  in  all  instances;  also  for  accurate 
measurements  and  for  truly  comparable  values,  the  correction 
should  be  applied  whenever  the  shift  exceeds  ±5  my.  However, 
to  produce  a  shift  in  the  cutoff  position  in  excess  of  5  my, 
a  temperature  difference  of  30  to  40  degrees  Celsius  is  neces¬ 
sary,  which  implies  that,  in  general,  a  correction  need  not 
be  applied  for  daily  variations  in  temperature,  and  only  in 
rare  cases  that  corrections  for  seasonal  temperature  varia¬ 
tions  are  necessary. 


. 
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The  solar  energy  above  the  upper  cutoff  wavelength  for 
the  filters  (which  is  at  about  2800  my  for  the  three  filters 
OG1,  RG2  and  RG8 )  reaching  the  earth's  surface  is  small  and  is 
a  function  of  the  precipitable  water  vapor  in  the  atmosphere. 
But  in  order  to  obtain  more  accurate  results  for  the  energy 
distribution  of  the  solar  radiation  reaching  the  earth's 
surface  within  different  spectral  bands,  a  small  correction 
AIL  should  be  added  to  the  filtered  radiation  being  reduced 
to  the  standard  conditions.  This  correction  is  small  and 
difficult  to  calculate,  so  it  is  often  omitted.  However,  the 
correction  AIL  has  been  calculated  as  a  function  of  the  pro¬ 
duct  of  air  mass  and  precipitable  water.  Table  3.1,  which 
is  a  reproduction  of  table  1,  Page  399,  Reference  [10],  lists 
this  correction. 


TABLE  3.1  Solar  intensity  above  the  longwave  cutoff  of  the 
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filters  AIL  (in  meal  .  cm  .  min  )  as  a  function  of 
atmospheric  humidity  nr  .w  (in  cm). 


m,  •  w 
h 

AIL 

m,  .  w 
h 

AIL 

m,  •  w 
h 

AIL 

0.10 

34 

0.60 

20 

4.30 

8 

0.13 

32 

0.80 

18 

6.30 

6 

0.16 

30 

1.00 

16 

11.00 

4 

0.21 

28 

1.10 

16 

22.00 

2 

0.26 

26 

1.50 

14 

44.0 

0 

0.34 

24 

2.10 

12 

0.45 

22 

3.00 

10 

\ 

■* 
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3.3  APPLICATION  OF  THE  FILTER  MEASUREMENTS 

The  method  of  dividing  the  solar  spectrum  into  two  regions 
by  means  of  each  filter  is  illustrated  by  considering  the 
red  filter  RG2 .  Its  wavelength  range  is  from  630  to  2800  my, 
and  F2  is  the  filter  factor.  If  the  intensity  of  the  total 
direct  solar  radiation  I  and  also,  practically  simultaneously, 
the  intensity  of  the  radiation  behind  the  filter  12  are  measured, 
then  by  multiplying  12  by  the  filter  factor  F2  and  adding 
AIL,  the  following  two  radiation  integrals  can  be  calculated 


oo 

I(X)dX  =  F2.I2  +  AIL  =  Ir  (3.1) 

630 

630 

» 

I ( X ) dX  =  I  -  Ir  =  Ik  (3.2) 

0 

where  I ( X )  represents  the  energy  distribution  within  the 
solar  spectrum, 

Ir  represents  the  red  and  infrared  radiation, 

Ik  represents  the  ultraviolet  and  visible  radiation. 

The  first  integral  (3.1)  covers  the  red  and  infrared  spectrum 
where  all  the  strong  absorption  bands  of  water  vapor  and  car¬ 
bon  dioxide  occur.  The  second  integral  (3.2)  covers  the  short 
wave  radiation,  where  there  is  almost  no  selective  absorption 
apart  from  the  ozone  absorption  in  high  layers  of  the  atmos¬ 
phere,  and  attenuation  of  the  radiation  in  this  region  is 


v 
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mainly  due  to  scattering  by  dry  gas,  dust  and  other  aerosols. 

If  II  and  18  are  the  intensities  measured  with  the  fil*- 
ters  0G1  and  RG8,  FI  and  F8  are  the  respective  filter  factors 
and  the  transmission  region  is  from  530  to  2800  my  for  OG1 
and  from  710  to  2800  my  for  RG3 ,  then  the  intensity  of  the 
solar  radiation  covered  by  the  following  integrals  can  be 
determined 


I(X)dX  =  FI- II  +  AIL  =  Id 

530 

530 


I ( X ) dX  =  I  -  Id  =  lb 

0 


(3.3) 


(3.4) 


and  I ( X) dX  =  F8-I8  +  AIL  =  Ii 

710 

710 


I (X) dX  =  I  -  Ii  =  Iv 

0 


(3.5) 


(3.6) 


where  Id 
lb 
Ii 
Iv 


represents  the  yellow  and  green  radiation, 
represents  the  blue  and  violet  radiation, 
represents  the  infrared  radiation  and 
represents  the  visible  radiation. 


In  general,  a  filter  is  used  to  separate  the  spectrum  of 
solar  radiation  into  two  well  defined  regions.  The  intensity 
of  the  radiation  below  the  short-wave  cutoff  of  the  filter  is 
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mainly  dependent  on  the  atmospheric  scattering  and  the  in¬ 
tensity  of  the  radiation  above  the  short-wave  cutoff  is  mainly 
dependent  on  the  selective  absorption  of  the  atmosphere. 

Data  concerning  the  variation  of  the  two  integrals  of  any  of 
the  three  filters  with  location  and  time  of  the  day  and  of 
the  year  is  of  value  not  only  for  studying  the  energy  distri¬ 
bution  in  specific  wavelength  ranges  but  also  for  studying 
the  constitution  of  the  atmosphere  with  regards  to  turbidity 
and  water  vapor  content.  The  latter  will  be  discussed  in 
chapter  four.  The  use  of  the  three  filters  0G1 ,  RG2  and  RG8 
allows  the  energy  distribution  to  be  calculated  for  seven  desig¬ 
nated  solar  spectrum  bands  as  shown  in  table  3.2. 


TABLE  3.2  Spectral  regions  of  the  solar  spectrum  obtained 


by  filter  measurements 


Denomenation 

Abbrevia¬ 

tion 

Wavelength  range 
(millimicrons) 

Calculation 

Blue  and  violet 

lb 

X  £  525 

I-Fl. 11 

Yellow  and  green 

Id 

525  <  X  <  630 

FI. I1-F2. 12 

Short-wave  (ultra¬ 
violet  and  visible) 

Ik 

X  £  630 

I-F2.I2 

Yellow  and  red 

ID 

525  <  X  <  710 

FI. I1-F8.I8 

Visible 

Iv 

X  <  710 

I-F8. 18 

Red  and  infrared 

Ir 

X  >  630 

F2. 12 

Infrared 

Ii 

X  >  710 

F8. 18 

. 


3.4  FILTER  MEASUREMENTS  AT  EDMONTON 
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The  output  of  the  three  normal  incidence  pyrheliometers 
mounted  on  a  tracking  unit  monitored  the  normal  incidence  insola¬ 
tion  over  the  observable  spectral  region  (i.e.  approximately 
0.3  to  3  microns)  and  the  insolation  passing  the  RG2  and  RG8 
filters  under  clear  sky  conditions  at  Edmonton.  This  allowed 
the  study  for  the  solar  energy  distribution  in  five  spectral 
bands.  The  procedures  of  preparing  the  hourly  records  of  the 
solar  energy  carried  within  the  total  normal  incidence  radia¬ 
tion  and  the  radiation  passing  RG2  and  RG8  filters  were 
described  in  section  six  of  chapter  one. 

The  filters  used  in  this  study  were  supplied  by  The  Eppley 
Laboratory.  They  were  2  millimeters  thick  and,  at  20  degrees 
Celsius,  had  a  center  of  lower  cutoff  at  628  my  for  RG2  filter 
and  692  my  for  RG8  filter  and  upper  cutoff  at  2800  my.  For 
calculations,  it  was  assumed  that  the  lower  center  of  cutoff 
was  at  626  my  for  RG2  filter  and  690  my  for  RG8  filter  which 
were  about  the  values  expected  at  the  average  yearly  outdoor 
temperature.  No  correction  was  applied  for  the  shift  in 
center  of  cutoff.  The  filter  factors  were  1.10  for  RG2  filter 
and  1.12  for  RG8  filter  and,  since  the  Eppley  normal  incidence 
pyrheliometers  have  a  quartz  window  between  the  radiation 
receiver  and  the  filter,  each  filter  factor  was  reduced  by 
one  percent  to  correct  for  multiple  reflection  between  the 
window  and  the  filter. 


The  correction  AIL  for  the  energy  above  2800  my  not 
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measured  by  the  filters  was  calculated  using  table  3.1.  This 
table  was  divided  into  two  intervals,  one  covering  the  range 
of  atmospheric  humidity  (precipitable  water  w  along  the  path  of 
the  sun's  rays)  for  which  m^*  w^  2  and  the  other  covering  the 
range  of  m^*  w  >  2.  Using  a  least  squares  technique  a  poly¬ 
nomial  was  fitted  for  each  interval.  The  correction  for  the 
solar  intensity  above  the  longwave  cutoff  of  the  filters  AIL 
was  calculated  using  the  polynomials: 

2  3 

AIL  =  33.61  -  33.25(mh*w)  +  19.76(mh*w)  -  0.22(mh*w)  , 

mhw  -  2  (3.7) 

2 

AIL  =  17.35  -  2.74(m^*w)  +  0.14(m^*w)  ,  m^w  >  2  (3.8) 

-2  -1 

where  AIL  is  m  meal. cm  .min  , 

(1  cal  cm  ^  min  =s  221.2  Btu  ft  ^hr  ^) 

w  is  the  precipitable  water,  and 

m,  is  the  relative  air  mass, 
h 

The  appropriate  relative  air  mass  values  were  calculated  and 
the  average  daily  precipitable  water  vapor  for  forenoon  and 
afternoon  were  calculated  using  the  published  values  in 
reference  [11] . 

The  use  of  the  two  filters  RG2  and  RG8  allowed  the  energy 
distribution  to  be  calculated  for  the  five  solar  spectral 
bands:  X  <  626,  X  >  626,  X  f  690,  X  >  690  and  626  <  X  1  690 
millimicrons.  Let  IDN,  12  and  18  be  the  measured  intensities 
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of  the  total  direct  normal  solar  radiation,  the  radiation 
passing  the  RG2  filter  and  the  radiation  passing  the  RG8  filter 
respectively  and  let  AIL  be  the  intensity  of  the  radiation 
above  the  long  wave  cutoff  of  the  filters.  The  intensities  of 
the  radiation  within  the  five  spectral  bands  were  calculated 
at  Edmonton  as  follows: 


Ir  =  (1.10  •  12)  +  AIL, 

Ik  =  IDN  -  Ir, 

Ii  =  (1.11  •  18)  +  AIL, 

Iv  =  IDN  -  Ii, 

Ie  =  Ir  -  Ii, 


X  >  62  6my 
X  -  626my 
X  >  6  9  0my 
X  <  6  90my 
626  <  X  1  690my 


where  Ir,  Ik,  Ii  and  Iv  are  as  defined  previously  and  Ie 
represents  the  red  band  radiation. 


3 . 5  THE  ANNUAL  VARIATION  OF  THE  SPECTRAL  DISTRIBUTION  OF  THE 
NORMAL  INCIDENCE  INSOLATION  MEASURED  AT  EDMONTON  AT 
DIFFERENT  HOUR  ANGLES 

Full  day,  half  day  and  part  day  data  were  used  in  the 
analysis.  For  each  forenoon  and  afternoon,  the  energy  within 
each  band  was  calculated  at  hour  angles  zero  to  six  when  it 
was  possible.  As  discussed  in  chapter  two,  the  normal  inci¬ 
dence  insolation  was  also  calculated  and  then  the  percentage 
of  the  normal  incidence  insolation  carried  by  the  radiation 


. 
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within  each  band  was  determined.  Figures  3.1  to  3.7  show  the 
latter  quantity  for  hour  angles  zero  to  six  throughout  the 
year.  The  curves  in  these  figures  were  best  fitted  by  eye  to 
show  the  approximate  annual  variation  of  the  fraction  of  the 
normal  insolation  with  the  five  spectral  bands.  The  fraction 
of  the  normal  insolation  for  the  three  bands  A  <  626my  (short¬ 
wave  fraction) ,  X  >  690my  (infrared  fraction)  and  626  <  X  < 

690my  (red  fraction)  should  be  read  on  the  left  hand  axis  from 
the  curves  named  Ikp,  lip  and  Iep  respectively,  and  the  frac¬ 
tions  for  the  two  bands  A  >  62  6my  and  A  <  6  90my  should  be 
read  on  the  right-hand  axis  from  the  curves  Ikp  and  lip, 
respectively. 

The  curves  in  figures  3.1  to  3.7  show  that  the  fraction 
lip  has  a  similar  annual  trend  at  different  hour  angles.  The 
same  type  of  results  are  also  found  for  the  fraction  Ikp.  The 
infrared  fractions  lip  have  flat  minima  during  the  summer 
months  June,  July  and  August  and  increase  gradually  during  the 
rest  of  the  year.  The  curves  representing  the  short  wave 
fractions  Ikp  have  flat  maxima  during  the  summer  months  and 
decrease  gradually  during  the  rest  of  the  year.  The  curves 
representing  the  red  fractions  Iep  show  slight  annual  varia¬ 
tion  at  all  hour  angles.  Table  3.3  lists  the  monthly  aver¬ 
age  of  the  infrared  fractions  and  the  shortwave  fractions  of 
the  normal  insolation.  The  fractions  of  the  normal  insolation 
within  the  two  bands  A  <  690my  (Ivp)  and  A  >  630my  (Irp)  can  be 
easily  determined  as  follows: 
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FIGURE  3.1  ANNUAL  VARIATION  OF  THE  FRACTIONS 

OF  IDN  WITHIN  DIFFERENT  WAVELENGTH 
•REGIONS  COMPUTED  AT  EDMONTON  AT 
HOUR  ANGLE  ZERO  . 


lip,  Ikp  and  Iep  (FRACTIONS  OF  IDN) 

.00  0.20  0.40v  0.60  0.80  /  1.00 
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FIGURE  3.2  ANNUAL  VARIATION  OF  THE  FRACTIONS 

OF  IDN  WITHIN  DIFFERENT  WAVELENGTH 
REGIONS  COMPUTED  AT  EDMONTON  AT 
HOUR  ANGLE  1  . 


.00  0.80  0.60,,  0.40  0.20  .  0.00 

irp  and  ivp  (FRACTIONS  OF  IDN) 


lip,  Ikp  and  Iep  (FRACTIONS  OF  I  DN ) 
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FIGURE  33  ANNUAL  VARIATION  OF  THE  FRACTIONS 

OF  I DN  WITHIN  DIFFERENT  WAVELENGTH 
.REGIONS  COMPUTED  AT  EDMONTON  AT 


HOUR  ANGLE  2  . 


.60.  0.40  0.20 

,p  (FRACTIONS  OF  IDN 


lip,  Ikp  and  Iep  (FRACTIONS  OF  IDN) 

.00  0.20  0.40'  0.60  0.80  '  1.00 
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FIGURE  3  -4  RNNURL  VARIATION  OF  THE  FRACTIONS 

OF  IDN  WITHIN  DIFFERENT  WAVELENGTH 


.REGIONS  COMPUTED  AT  EDMONTON  AT 


HOUR  ANGLE  3. 


.  60  0.40  0.20  . 

,p  (FRACTIONS  OF  I DN) 
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FIGURE  3.5  ANNUAL  VARIATION  OF  THE  FRACTIONS 

OF  IDN  WITHIN  DIFFERENT  WAVELENGTH 
•REGIONS  COMPUTED  AT  EDMONTON  AT 
HOUR  ANGLE  4-  , 
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FIGURE  3.6  ANNUAL  VARIATION  OF  THE  FRACTIONS 

OF  IDN  WITHIN  DIFFERENT  WAVELENGTH 
REGIONS  COMPUTED  AT  EDMONTON  AT 
HOUR  ANGLE  5. 
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FIGURE  3-7 


ANNUAL  VARIATION  OF  THE  FRACTIONS 
OF  ION  WITHIN  DIFFERENT  WAVELENGTH 
REGIONS  COMPUTED  AT  EDMONTON  AT 


HOUR  ANGLE  6. 


).  0.40  0.20  . 

(FRACTIONS  OF  idn) 
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TABLE  3  »3  Monthly  average  percentage  of  the  normal  incidence 

radiation  in  the  infrared  spectrum  band  lip ( A>690my) 


and  the  shortwave  spectrum  band  Ikp(A<626my)  at 


Edmonton  at  different  hour  angles. 


Month 

Hour  Angle  From  Solar  Noon 

0 

1 

2 

3 

4 

5 

6 

January 

lip 

64.6 

64.6 

66.2 

73.0 

Ikp 

28.7 

28.9 

27.6 

21.5 

February 

lip 

61.2 

61.3 

63.0 

66.3 

70.6 

Ikp 

32.2 

32.1 

30.2 

26.7 

21.0 

March 

lip 

59.2 

59.3 

59.7 

61.7 

65.3 

68.7 

Ikp 

34.8 

34.5 

33.8 

31.5 

28.1 

26.1 

April 

lip 

56.4 

56.6 

57.8 

58.8 

61.2 

65.4 

67.4 

Ikp 

36.0 

35.9 

35.2 

34.6 

32.8 

28.4 

24.6 

May 

lip 

55.6 

55.7 

55.8 

56.1 

57.4 

59.6 

65.6 

Ikp 

37.6 

37.2 

36.8 

36.4 

34.2 

31.0 

26.4 

June 

lip 

54.1 

54.3 

54.5 

55.2 

55.9 

58.1 

62.8 

Ikp 

37.5 

37.2 

37.1 

36.3 

35.7 

32.9 

30.2 

July 

lip 

54.1 

54.5 

54.6 

55.5 

56.4 

58.3 

62.2 

Ikp 

39.9 

39.9 

39.7 

39.4 

38.4 

35.8 

31.5 

August 

lip 

54.1 

54.5 

54.9 

55.3 

56.0 

59.9 

67.2 

Ikp 

39.0 

38.2 

38.0 

37.4 

35.2 

32.8 

25.2 

September 

lip 

55.0 

55.2 

55.4 

56.3 

58.1 

62.9 

69.7 

Ikp 

38.6 

38.4 

38.2 

38.2 

35.9 

31.8 

24.2 

October 

lip 

58.7 

58.9 

59.6 

62.1 

64.6 

Ikp 

35.3 

35.2 

34.5 

31.8 

29.0 

November 

lip 

61.5 

61.6 

63.2 

65.7 

70.3 

Ikp 

33.4 

32.9 

31.7 

27.8 

23.7 

December 

lip 

65.5 

66.2 

70.3 

73.7 

Ikp 

28.8 

27.1 

24.2 

20.2 

v 
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Ivp  =  100  -  lip 
and  Irp  =  100  -  Ikp 

Figure  3.8  depicts  all  the  curves  in  figures  3.1  to  3.7 
on  one  graph  for  visual  comparison  at  different  hour  angles. 
The  infrared  fraction  curves  have  close  agreement  at  hour 
angles  0,  1  and  2  with  an  average  value  of  54  percent  in  May, 
June,  July  and  August  with  a  variation  of  ±1  percent  about  the 
average.  This  average  increases  to  about  62  percent  for  the 
rest  of  the  year  with  a  variation  of  ±6  percent  about  the 
average.  As  the  hour  angle  exceeds  2  the  infrared  fractions 
increase,  a  characteristic  which  occurs  for  all  days  through¬ 
out  the  year.  Close  agreement  between  the  curves  up  to  hour 
angle  4  is  noticed  during  the  summer  months,  an  average  value 
of  about  55  percent  can  be  assumed  for  May  to  August  at  hour 
angles  0  to  4 .  At  hour  angle  5  and  6  the  solar  energy  is 
available  for  a  short  period  of  the  year,  approximately  from 
April  to  September,  and  the  infrared  fractions  have  an  average 
value  of  about  66  percent. 

The  shortwave  fraction  curves  show  that  the  maximum 
fraction  of  about  38  percent  occurs  in  the  summer  at  hour 
angles  0,  1  and  2  then  the  fractions  decrease  on  both  sides 
of  the  curve  to  about  31±6  percent  for  the  rest  of  the  year. 

As  the  hour  angle  exceeds  2  the  shortwave  fractions  decrease, 
a  characteristic  which  occurs  for  all  days  throughout  the 
year.  At  hour  angle  six  the  shortwave  fractions  reach  a  value 
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FIGURE  3.3  ANNUAL  VARIATION  OF  THE  FRACTIONS 

OF  IDN  WITHIN  DIFFERENT  WAVELENGTH 
REGIONS  COMPUTED  AT  EDMONTON  AT 


HOUR  ANGLES  O  to  6 . 


(FRACTIONS  OF  IDN) 
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of  25±5  percent. 

The  fractions  within  the  red  band  have  slight  variation 
at  different  hour  angles  for  all  days  throughout  the  year  and 
the  average  value  is  about  6.5  percent. 

The  variation  of  the  fractions  of  the  normal  insolation 
in  different  bands  with  time  of  the  day  and  the  year  shows  an 
agreement  with  the  fact  that  the  attenuation  of  the  solar  beam 
due  to  dry  air  molecules,  water  vapor,  dust  and  other  aerosols 
in  the  atmosphere  is  not  uniform  over  the  solar  spectrum. 

The  closer  agreement  between  the  infrared  fraction  curves 
at  different  hour  angles  till  hour  angle  4  during  the  summer 
months  is  due  to  relatively  small  change  in  the  path  length 
of  the  radiation  through  the  atmosphere  while  the  precipi- 
table  water,  which  is  the  major  attenuating  factor  of  the  in¬ 
frared  radiation,  is  at  about  the  same  level  during  the  summer. 

The  decrease  of  the  shortwave  fractions  as  the  hour  angle 
increases  and  also  in  the  winter  months  is  due  to  the  longer 
path  through  the  atmosphere  where  more  attenuation  is  caused 
by  scattering  by  dust  and  other  aerosols,  which  are  the  major 
attenuating  factors  of  the  shortwave  radiation,  as  well  as 
scattering  by  air  molecules  which  attenuate  the  shortwave 
radiation  more  than  the  infrared  radiation. 

The  monthly  averages  at  different  hour  angles  of  the  mag¬ 
nitude  of  the  energy  within  each  of  the  five  spectral  bands 
calculated  at  Edmonton  are  listed  in  tables  3.4  t.o  3.6. 
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TABLE  3 . 4  Monthly  average  normal  incidence  solar  radiation 
in  the  red  and  infrared  band  Ir  (A>626my)  and  the 

shortwave  band  Ik  (A<626my)  measured  at  Edmonton  at 

-2  -1 

different  hour  angles.  Unit:  Btu.ft  .hr 


Month 

Hour 

Angle 

From  Solar  Noon 

0 

1 

2 

3 

4 

5 

6 

January 

Ir 

181 

178 

164 

139 

Ik 

73 

72 

63 

40 

February 

Ir 

192 

188 

183 

170 

149 

1 

Ik 

91 

89 

79 

62 

41 

March 

Ir 

190 

189 

181 

178 

158 

138 

Ik 

101 

99 

94 

82 

64 

45 

April 

Ir 

189 

186 

183 

179 

171 

150 

127 

Ik 

106 

105 

102 

94 

83 

60 

29 

May 

Ir 

185 

182 

181 

176 

175 

159 

133 

Ik 

104 

110 

105 

102 

90 

75 

48 

June 

Ir 

184 

182 

182 

178 

175 

159 

145 

Ik 

110 

108 

107 

102 

97 

79 

62 

July 

Ir 

177 

174 

172 

170 

162 

153 

143 

Ik 

113 

112 

110 

100 

94 

86 

66 

August 

Ir 

174 

172 

170 

164 

155 

141 

115 

Ik 

120 

118 

104 

105 

92 

73 

48 

September 

Ir 

178 

178 

171 

168 

155 

132 

111 

Ik 

112 

111 

103 

106 

95 

65 

35 

October 

Ir 

179 

178 

174 

165 

151 

126 

Ik 

106 

104 

99 

84 

64 

42 

November 

Ir 

174 

171 

163 

146 

129 

Ik 

89 

85 

76 

57 

40 

December 

Ir 

171 

168 

153 

129 

Ik 

66 

62 

48 

32 

.  - 1 
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TABLE  3 .5  Monthly  average  normal  incidence  solar  radiation 
in  the  infrared  band  Ii  (A>690my)  and  the  visible 

band  Iv  (A<690my)  measured  at  Edmonton  at  different 

-2  -1 

hour  angles.  Unit:  Btu.ft  .hr 


Month 

Hour 

Angle 

From  Solar  Noon 

0 

1 

2 

3 

4 

5 

6 

January 

Ii 

164 

161 

150 

121 

Iv 

90 

89 

77 

63 

February 

Ii 

168 

166 

162 

153 

135 

Iv 

109 

107 

97 

79 

55 

March 

Ii 

173 

169 

163 

161 

143 

128 

Iv 

119 

116 

112 

100 

77 

58 

April 

Ii 

165 

167 

165 

160 

154 

138 

114 

Iv 

131 

130 

121 

113 

98 

73 

43 

May 

Ii 

162 

160 

158 

154 

148 

136 

111 

Iv 

130 

129 

126 

122 

111 

94 

67 

June 

Ii 

158 

155 

154 

152 

147 

137 

126 

Iv 

137 

133 

130 

126 

118 

100 

81 

July 

Ii 

158 

156 

153 

152 

147 

132 

130 

Iv 

132 

131 

129 

126 

115 

99 

79 

August 

Ii 

155 

153 

151 

149 

139 

128 

107 

Iv 

132 

131 

126 

122 

109 

89 

54 

September 

Ii 

161 

160 

159 

155 

137 

122 

101 

Iv 

130 

127 

121 

116 

105 

78 

44 

October 

Ii 

165 

164 

160 

153 

130 

110 

Iv 

121 

115 

112 

98 

80 

43 

November 

Ii 

157 

155 

151 

135 

119 

Iv 

99 

95 

89 

64 

50 

December 

Ii 

154 

152 

142 

120 

Iv 

82 

78 

59 

42 

. 
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TABLE  3.6  Monthly  average  normal  incidence  solar  radiation 

in  the  wavelength  range  of  626  to  690  millimicrons 

(red  band  Ie)  measured  at  Edmonton  at  different 

-2  -1 

hour  angles.  Unit:  Btu-ft  *hr 


Month 

Hour 

Angle  From  Solar  Noon 

0 

1 

2 

3 

4 

5 

6 

January 

17 

16 

14 

12 

February 

19 

19 

18 

14 

11 

March 

19 

19 

18 

14 

12 

7 

April 

18 

16 

15 

12 

11 

9 

8 

May 

16 

14 

14 

13 

12 

11 

10 

June 

16 

15 

15 

14 

12 

11 

10 

July 

18 

17 

16 

15 

14 

13 

12 

August 

17 

17 

16 

16 

15 

14 

12 

September 

18 

18 

18 

17 

15 

12 

9 

October 

18 

17 

17 

16 

15 

13 

November 

16 

14 

12 

10 

8 

December 

14 

14 

12 

9 

«* 
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The  maximum  energy  within  the  red  and  infrared  band  Ir 
and  the  infrared  band  Ii  occurs  in  March  and  April  when  the 
precipitable  water  vapor  in  the  atmosphere  is  low.  Lower 
energy  is  available  in  July  and  August  when  the  atmosphere  has 
a  high  precipitable  water  level.  The  least  energy  occurs 
during  the  winter  months  at  large  hour  angles  (when  the  solar 
altitude  is  less  than  20  degrees)  . 

For  the  shortwave  band  Ik  and  the  visible  band  Iv  the 
maximum  energy  occurs  in  the  summer  months  May  to  August  at 
hour  angles  0  to  3  when  the  path  length  of  the  radiation  through 
the  atmosphere  is  relatively  short.  However,  the  magnitude 
of  the  energy  is  slightly  larger  than  the  values  measured  during 
April  and  September  when  the  path  length  through  the  atmos¬ 
phere  is  relatively  longer.  This  is  probably  due  to  an  increase 
of  the  dust  particles  in  the  atmosphere  during  May  to  August. 

The  least  energy  is  measured  during  the  winter  months  at 
large  hour  angles.  The  magnitude  of  the  energy  within  the  red 
band  Ie  has  a  similar  trend. 


*■ 
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3.6  VARIATION  OF  THE  INFRARED  AND  SHORTWAVE  FRACTIONS  OF  THE 

NORMAL  INSOLATION  WITH  THE  SOLAR  ALTITUDE 

The  variation  of  the  infrared  fraction  and  the  shortwave 
fraction  of  the  direct  normal  incidence  solar  radiation  (normal 
insolation)  was  studied.  The  average  fraction  of  the  energy 
within  each  of  the  infrared  band  and  the  shortwave  band  was 
calculated  for  each  one  degree  interval  of  change  in  solar 
altitude  from  5  to  60  degrees.  Figures  3.9  and  3.10  show  the 
plot  of  these  values.  It  seems  that  the  infrared  fraction 
increases  somewhat  uniformly  from  about  54  percent  at  60  degrees 
to  about  62  percent  at  20  degrees  and  then  increases  rapidly 
to  about  74  percent  at  5  degrees  solar  altitude. 

Also,  it  appears  that  the  shortwave  fraction  decreases 
somewhat  uniformly  from  about  40  percent  at  60  degrees  to  about 
35  percent  at  20  degrees  and  then  decreases  to  about  20  percent 
at  5  degrees  solar  altitude. 
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Figure  3.9  Variation  of  the  infrared  fraction,  lip,  of  the  direct 
normal  incidence  solar  radiation,  IDN,  with  the  solar 
altitude  (degrees). 
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Figure  3.10  Variation  of  the  shortwave  fraction,  Ikp,  of  the  direct 
normal  incidence  solar  radiation,  IDN,  with  the  solar 
altitude  (degrees). 
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3.7  VARIATION  OF  THE  INFRARED  AND  SHORTWAVE  FRACTIONS  OF  THE 

NORMAL  INSOLATION  WITH  THE  OPTICAL  AIR  MASS 

The  attenuation  of  solar  radiation  is  a  function  of  the 
path  length  through  the  atmosphere  and  each  station,  due  to 
elevation,  has  a  different  path  length.  To  compensate  for  the 
effect  of  elevation,  calculations  are  performed  using  the 
optical  air  mass.  The  optical  air  mass  is  obtained  by  multi¬ 
plying  the  relative  air  mass  (which  is  approximately  equal  to 
1/sin  (altitude))  by  the  ratio  of  the  station  pressure  to 
standard  atmospheric  pressure.  Representation  of  the  frac¬ 
tions  of  the  normal  insolation  as  functions  of  optical  air 
mass  seems  preferable. 

From  the  results  of  the  annual  variation  of  the  infrared 
and  shortwave  fractions  and  their  variation  with  solar  alti¬ 
tude,  it  was  decided  to  represent  the  fractions  as  functions 
of  the  optical  air  mass  at  Edmonton  for  two  periods  of  the 
year  and  two  intervals  of  solar  altitude.  The  first  period 
includes  April  to  September  and  the  second  period  includes 
October  to  March.  The  solar  altitude  intervals  are  from  5 
to  20  degrees  and  from  20  to  60  degrees.  A  study  of  the 
variation  of  the  solar  altitude  with  time  of  the  day  and 
time  of  the  year  at  Edmonton  showed  that  for  the  first  period, 
the  solar  altitude  is  greater  than  20  degrees  on  all  days 
till  hour  angle  5  except  in  April,  the  second  half  of  August 
and  September  where  it  is  valid  only  till  hour  angle  4.  The 
solar  altitude  is  less  than  20  degrees  after  hour  angle  5 


- 

* 
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through  out  the  first  period.  For  the  second  period,  the 
solar  altitude  is  less  than  20  degrees  except  in  March  from 
hour  angle  zero  to  three  and  in  February  and  October  till 
hour  angle  2. 

In  Figures  3.11  to  3.14  the  shortwave  and  the  infrared 
fractions  are  represented  as  functions  of  the  optical  air  mass 
for  the  two  solar  altitude  intervals  5  to  20  degrees  and  20 
to  60  degrees  for  both  periods  April  to  September  and  October 
to  March.  The  average  optical  air  mass  was  calculated  for 
each  one  degree  interval  of  solar  altitude.  It  appeared  that 
the  fractions  could  be  represented  by  linear  relationships 
and  the  best  fit  line  was  calculated  for  each  set  of  data.  The 
resulting  equations  representing  the  infrared  fractions  as 
functions  of  the  optical  air  mass  m  are  listed  in  table  3.7  for 
the  two  intervals  of  solar  altitude  and  the  two  periods  of 
the  year.  For  solar  altitude  angles  a  >  20°  the  relative  air 
mass  is  given  by  l/(sin  a)  and  assuming  that  the  optical  air 
mass  is  calculated  using  mean  annual  station  pressure  (933 
millibars) ,  approximate  equations  can  be  written  in  terms  of 
the  solar  altitude  a  for  a  >  20°.  These  equations  are  also 
listed  in  table  3.7.  Similar  linear  relationships  for  the 
shortwave  fractions  are  given  in  table  3.8. 

For  the  same  optical  air  mass,  the  infrared  fraction  is 
less  in  the  summer  than  in  the  winter  by  about  2.5  percent. 

On  the  other  hand,  the  shortwave  fraction  is  about  1.5  percent 
higher  in  the  summer. 


* 
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OPTICAL  AIR  MASS 


Figure  3.11  Variation  of  the  infrared  fraction,  lip,  of  the  direct 
normal  incidence  solar  radiation,  IDN,  with  the  optical 
air  mass  for  the  solar  altitude  interval  20  to  60  degrees 
throughout  the  year. 
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OPTICAL  AIR  MASS 


Figure  3.12  Variation  of  the  infrared  fraction,  lip,  of  the  direct 
normal  incidence  solar  radiation,  IDN,  with  the  optical 
air  mass  for  the  solar  altitude  interval  5  to  20  degrees 
throughout  the  year. 
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Figure  3.13  Variation  of  the  shortwave  fraction,  Ikp,  of  the  direct 
normal  incidence  solar  radiation,  IDN,  with  the  optical 
air  mass  for  the  solar  altitude  interval  20  to  60  degrees 
throughout  the  year. 
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Figure  3.14  Variation  of  the  shortwave  fraction,  Ikp,  of  the  normal 

incidence  solar  radiation,  IDN,  with  the  optical  air  mass 
for  the  solar  altitude  interval  5  to  20  degrees  throughout 
the  year. 
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TABLE  3 . 7  Equations  representing  the  infrared  fraction 

(A  >  690my)  of  the  normal  incidence  solar 
radiation  in  terms  of  the  optical  air  mass  m 
and  the  solar  altitude  a. 


(*) 

Equation 

Altitude 

Interval 

Time  of  the  Year 

Fraction  =  0.513  +  0.030  •  m 

20  <  a  <  60 

April  to  September 

=  0.513  +  0.028  •  (1/sina) 

II 

II 

Fraction  =  0.528  +  0.033  •  m 

II 

October  to  March 

=  0.528  +  0.031  •  (1/sina) 

II 

II 

Fraction  =  0.555  +  0.019  •  m 

5  <  a  <  20 

April  to  September 

Fraction  =  0.585  +  0.017  •  m 

II 

October  to  March 

(*) 

The  straight  lines  are  shown  in  figures  3.11  and  3.12. 
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TABLE  3 . 8  Equations  representing  the  shortwave  fraction 

(A  <  626my)  of  the  normal  incidence  solar  radiation 
in  terms  of  the  optical  air  mass  m  and  the  solar 
altitude  a. 


(*) 

Equation  ' 

Altitude 

Interval 

Time  of  the  Year 

Fraction  =  0.420  -  0.031  •  m 

20  <  a  <  60 

April  to  September 

=  0.420  -  0.029  •  (1/sina) 

II 

II 

Fraction  =  0.451  -  0.051  •  m 

II 

October  to  March 

=  0.451  -  0.048  •  (1/sina) 

II 

II 

Fraction  =  0.373  -  0.017  •  m 

5  <  a  <  20 

April  to  September 

Traction  =  0.352  -  0.017  •  m 

II 

October  to  March 

(*> 

The  straight  lines  are  shown  in  figures  3.13  and  3.14. 
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3 . 8  DETERMINATION  OF  ATMOSPHERIC  TRANSMISSION  FACTORS  AND 

EXTINCTION  COEFFICIENTS  FOR  DIRECT  RADIATION  WITHIN  THE 

WHOLE  SOLAR  SPECTRUM  AND  THE  FIVE  SPECTRAL  BANDS 

3.8.1  THE  ATMOSPHERIC  TRANSMISSION  FACTORS:  As  mentioned 
in  section  2.3  the  atmospheric  transmission  factor  for  the 
integral  wavelength  direct  normal  incidence  radiation  can  be 
calculated,  using  a  simplified  formula  of  equation  2.6,  by 
dividing  the  measured  radiation  at  the  earth's  surface  reduced 
to  mean  earth  -  sun  distance  by  the  solar  constant.  The 
atmospheric  transmission  factors  for  a  given  band  of  the  spec¬ 
trum  can  be  calculated  by  dividing  the  measured  radiation  within 
this  band  reduced  to  mean  earth-sun  distance  by  the  portion  of 
the  solar  constant  within  the  same  band. 

The  portion  of  the  solar  constant  within  each  of  the 
five  spectral  bands  under  investigation  in  this  study  was 
calculated  from  table  1  in  reference  [33] .  The  intensities  of 
the  direct  radiation  within  the  whole  spectrum  (IDN) ,  infrared 
(Ii) ,  visible  (Iv) ,  red  and  infrared  (Ir) ,  shortwave  (Ik)  and 
red  (Ie)  bands  were  reduced  to  mean  earth-sun  distance.  The 
atmospheric  transmission  factors  TDN,  Ti ,  Tv,  Tr,  Tk  and  Te 
for  the  whole  spectrum  and  the  five  bands  were  then  calculated. 

Tables  3.9  and  3.10  list  the  average  values  for  the 
energy  quantities  IDN,  Ii,  Iv,  Ir,  Ik  and  le  reduced  to  mean 
earth-sun  distance,  the  average  optical  air  mass  m  and  the 
average  values  for  the  transmission  factors  TDN,  Ti,  Tv,  Tr, 
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TABLE  3 . 9  Normal  incidence  solar  radiation,  IDN ,  at  different 
solar  altitudes  a  (degrees)  and  its  distribution  in 


different  spectral  bands  throughout  the  year. 
Uni±:  Btu.ft  ^.hr  ^ . 


— 

a 

(*) 

Period 

m 

-2  -1 

Normal  incidence  radiation  (Btu.ft  .hr  ) 

IDN 

Ii 

V 

Iv 

Ir 

Ik 

Ie 

60 

1 

1.09 

296 

160 

136 

174 

118 

18 

50 

1 

1.21 

290 

159 

131 

177 

113 

17 

40 

1 

1.48 

289 

162 

127 

179 

109 

17 

2 

1.45 

289 

163 

126 

181 

108 

17 

30 

1 

1.88 

264 

150 

113 

168 

96 

17 

2 

1.85 

283 

169 

114 

181 

101 

12 

20 

1 

2.77 

229 

137 

91 

153 

76 

15 

2 

2.72 

262 

162 

99 

178 

83 

15 

15 

1 

3.63 

210 

132 

77 

146 

63 

14 

2 

3.62 

243 

158 

84 

171 

72 

12 

10 

1 

5.32 

157 

105 

52 

117 

39 

12 

2 

5.70 

203 

140 

62 

152 

51 

8 

(*) 


Period  1  includes  the  months  April  to  September. 
Period  2  includes  the  months  October  to  March. 
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TABLE  3.10  Atmospheric  transmission  factors  for  the  whole  solar 
spectrum  and  the  five  bands  calculated  at  different 
solar  altitudes  a  (degrees)  throughout  the  year. 


a 

(*) 

Period 

m 

Atmospheric 

transmission 

factors 

TDN 

Ti 

Tv 

Tr 

Tk 

Te 

60 

1 

1.09 

0.690 

0.717 

0.691 

0.702 

0.711 

0.591 

50 

1 

1.21 

0.676 

0.715 

0.663 

0.699 

0.679 

0.579 

40 

1 

1.48 

0.674 

0.725 

0.646 

0.708 

0.660 

0.577 

2 

1.45 

0.674 

0.731 

0.640 

0.714 

0.651 

0.577 

30 

1 

1.88 

0.615 

0.675 

0.577 

0.663 

0.577 

0.580 

2 

1.85 

0.660 

0.756 

0.581 

0.716 

0.612 

0.415 

20 

1 

2.77 

0.534 

0.616 

0.565 

0.603 

0.458 

0.511 

2 

2.72 

0.611 

0.729 

0.504 

0.704 

0.504 

0.515 

15 

1 

3.63 

0.489 

0.593 

0.395 

0.578 

0.382 

0.469 

2 

3.62 

0.566 

0.711 

0.428 

0.675 

0.434 

0.405 

10 

1 

5.32 

0.366 

0.470 

0.265 

0.463 

0.239 

0.413 

2 

5.70 

0.473 

0.631 

0.316 

0.602 

0.305 

0.291 

{ *) 

'  Period  1  includes  the  months  April  to  September. 
Period  2  includes  the  months  October  to  March. 
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Tk  and  Te  for  the  two  periods  of  the  year,  April  to  September 
and  October  to  March,  at  solar  altitudes  60,  50,  40,  30,  20, 

15  and  10  degrees.  In  general,  for  the  same  solar  altitude, 
the  atmospheric  transmission  factors  are  lower  for  the  period 
April  to  September  when  the  path  length  through  the  atmosphere, 
optical  air  mass,  has  higher  level  of  precipitable  water  and 
more  dust  particles. 

3.8.2  THE  EXTINCTION  COEFFICIENTS:  As  mentioned  in 
section  2.4  Bouguer's  law  (equation  2.1)  may  be  used  to  represent 
the  attenuation  of  the  normal  incidence  radiation  within  the 
whole  solar  spectrum.  Also,  it  may  be  used  to  represent  the 
normal  incidence  radiation  within  different  spectral  bands  [27]. 
The  portion  of  the  solar  constant  within  the  band  of  concern 
should  be  used  in  applying  Bouguer's  law  to  that  band. 

Bouguer's  law  was  applied  to  the  measured  energy  quantities 
IDN,  Ii,  Iv,  Ir,  Ik  and  Ie  reduced  to  mean  earth-sun  distance 
and  the  extinction  coefficients  for  the  total  spectrum  (ADN) , 
the  infrared  band  (Ai) ,  visible  band  (Av) ,  red  and  infrared 
band  (Ar) ,  shortwave  band  (Ak)  and  the  red  band  (Ae)  were 
calculated  at  different  hour  angles  throughout  the  year.  The 
average  value  for  each  of  these  coefficients  were  calculated 
for  the  two  periods  of  the  year,  April  to  September  and  October 
to  March,  at  solar  altitudes  60,  50,  40,  30,  20,  15  and  10 
degrees.  Table  3.11  lists  these  values  and  the  corresponding 
average  optical  air  mass  m. 


f. 
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TABLE  3.11  Extinction  coefficients  for  the  whole  solar  spec¬ 
trum  and  the  five  bands  calculated  at  different 
solar  altitudes  a  (degrees)  throughout  the  year. 


a 

( * ) 

Period  v  ' 

m 

Extinction 

Coefficients 

ADN 

Ai 

Av 

Ar 

Ak 

Ae 

60 

1 

1.09 

0.340 

0.306 

0.338 

0.325 

0.312 

0.482 

50 

1 

1.21 

0.234 

0.277 

0.340 

0.296 

0.319 

0.451 

40 

1 

1.48 

0.267 

0.217 

0.295 

0.233 

0.281 

0.371 

2 

1.45 

0.272 

0.216 

0.308 

0.233 

0.296 

0.379 

30 

1 

1.88 

0.258 

0.209 

0.293 

0.218 

0.292 

0.290 

2 

1.85 

0.225 

0.151 

0.294 

0.181 

0.265 

0.476 

20 

1 

2.77 

0.277 

0.175 

0.276 

0.182 

0.282 

0.242 

2 

2.72 

0.181 

0.116 

0.252 

0.129 

0.252 

0.244 

15 

1 

3.63 

0.197 

0.144 

0.256 

0.151 

0.265 

0.209 

2 

3.62 

0.157 

0.094 

0.234 

0.109 

0.231 

0.250 

10 

1 

5.32 

0.189 

0.142 

0.250 

0.145 

0.269 

0.166 

2 

5.70 

0.131 

0.081 

0.202 

0.089 

0.208 

0.217 

^ Period  1  includes  the  months  April  to  September. 
Period  2  includes  the  months  October  to  March. 


91 


In  general,  at  the  same  solar  altitude,  the  extinction 
coefficients  are  higher  for  the  period  April  to  September 
than  for  the  period  October  to  March.  The  ratio  is  1.29  for 
the  whole  spectrum,  1.47  for  the  infrared  band  and  1.16  for 
the  shortwave  band.  Also,  it  was  found  that: 


(Ak/ADN)  =  1.15 

(Ak/ADN) 2 

=  1.4 

(Ai/ADN) 1  =  0.79 

(Ai/ADN) 2 

=  0.66 

(Ak/Ai) x  =  1.45 

(Ak/Ai) 2 

=  2.1 

where  subscript  1  refers  to  the  period  April  to  September  and 
subscript  2  refers  to  the  period  October  to  March. 


‘ 

CHAPTER  4 


TURBIDITY  OF  THE  ATMOSPHERE  AND 
ABSORPTION  OF  SOLAR  RADIATION 
BY  ATMOSPHERIC  WATER  VAPOR 


92 


93 


4.1  ATMOSPHERIC  TURBIDITY 

The  basic  relation  for  the  depletion  of  solar  radiation 
by  the  earth's  atmosphere  was  given  in  chapter  2  by: 

GO 

2  -A ( X ) *m 

I  =  (RQ/Ri)  1 0  ( A )  •  e  •  dX  (4.1) 

o 

Rewriting  equation  4.1  in  terms  of  transmission  of  the  dif¬ 
ferent  components  of  the  earth's  atmosphere  gives: 


• 

H 

II 

00 

* 

I0(A)  •  tRm  (A)  •  xDm(A)  •  t”(A)  •  x0m(A)  • 

o 

tCq  ( A )  ■  dA  (4.2) 

where  IQ(X) 

is  the  solar  constant  as  a  function  of  the 

wavelength  X, 

I 

is  the  measured  solar  intensity  at  the  ground. 

R 

=  (R-l/R0)2,  the  reduction  factor  for  mean 

earth-sun  distance, 

R. 

1 

is  earth-sun  distance  at  time  of  observation, 

R 

o 

mean  earth-sun  distance. 

trU> 

transmission  of  a  Rayleigh  atmosphere,  clean 

and  dry  atmosphere,  of  unit  air  mass, 

tdU) 

transmission  due  to  dust  in  a  unit  air  mass, 

VA) 

transmission  due  to  water  vapor  in  a  unit  air 

mass , 

T0U) 

transmission  due  to  ozone  in  a  unit  air  mass, 
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t  (A)  transmission  due  to  carbon  dioxide  in  a  unit 
air  mass,  and 

m  is  the  optical  air  mass. 

The  turbidity  of  the  atmosphere  is  defined  as  the  redu¬ 
ced  transmission  of  the  atmosphere,  caused  by  absorption  and 
scattering  by  solid  or  liquid  particles,  other  than  clouds, 
held  in  suspension.  On  the  basis  of  the  physical  considera¬ 
tions  several  methods  for  measuring  atmospheric  turbidity 
have  been  proposed.  Angstrom  (1929,  1930)  defined  a  turbidity 
coefficient  3/  the  extinction  coefficient  at  A  =  ly  and  y, 
a  wavelength  exponent.  Another  turbidity  coefficient  B,  was 
defined  by  Schiiepp  (1949)  to  represent  the  extinction  at  A  = 
0.5y.  A  third  measure  of  turbidity,  called  Linke's  turbidity 
factor  T,  was  proposed  by  Linke  as  early  as  1922. 

4.1.1  THE  TURBIDITY  FACTOR,  T:  It  is  defined  as  the 
ratio  of  the  extinction  of  direct  solar  radiation  in  the  real 
atmosphere  to  that  of  a  Rayleigh  atmosphere,  i.e.,  an  atmos¬ 
phere  free  of  dust,  water  vapor  and  ozone. 

In  1922  Linke  derived  a  mathematical  formula  for  a 
separation  of  the  extinction  by  dust  and  the  selective  absorp¬ 
tion  by  the  gases  of  the  atmosphere  from  the  scattering  by 
the  molecules.  He  proposed  that: 


1^  =  IQX  •  exp(-aR(A)  •  T  •  m) 


(4.3) 
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where  a  (A) 
R 


i 


OA 


I 


A 


is  the  extinction  coefficient  of  molecular 
scattering  for  clear  dry  atmosphere, 
is  the  intensity  of  parallel  monochromatic 
beam  at  air  mass  zero, 

is  the  intensity  of  parallel  monochromatic 
beam  at  the  surface. 


Linke's  turbidity  factor  T  expresses  the  total  extinction  with 
the  molecular  scattering  as  unity.  For  an  ideal  pure  (clean 
and  dry)  atmosphere,  T  becomes  one.  When  T  equals  2.0,  the 
extinction  of  the  real  atmosphere  containing  water  vapor  and 
dust  is  twice  that  of  a  perfectly  clean  atmosphere  with  no 
selective  absorption.  Therefore  the  value  of  T  can  never  be 
smaller  than  one. 

By  integrating  equation  4.3  and  applying  the  reduction 
for  mean  solar  distance,  the  total  radiation  of  all  wave¬ 
lengths  is  given  by: 


I 


1/R  • 


exp (-a  (A)  *  T  •  m)  •  dA 
R 


R  • I  =  Iq  •  exp(-T  •  a  (m)  •  m) 


(4.5) 


where  Iq 


oo 

Iq(A)  dA  represents  the  solar  constant, 

0 


a  (m)  is  the  mean  value  over  all  wavelengths  of  the 
R 

extinction  coefficient  in  clean  dry  atmosphere, 
and  is  now  a  function  of  air  mass  m,  and 
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T  is  the  turbidity  factor  for  total  radiation  of  all 
wavelengths . 

From  equation  4.5  a  formula  for  the  computation  of  the  tur¬ 
bidity  factor  is  deduced  as: 

T  =  P (m)  •  [log(IQ)  -  log(I)  -  log(R)]  (4.6) 

where  P (m)  =  (m  •  a  (m)  •  log(e))  \  and  is  given  in  table  6  of 

K. 

the  I GY  Instruction  Manual  1958  [10]. 

The  optical  air  mass  is  used  in  evaluation  of  the  turbi¬ 
dity  factor  T.  For  high  level  stations,  when  tables  are 
used  for  the  computation  of  turbidity,  the  air  mass  m  =  1  is 
not  the  vertical  air  mass,  with  the  sun  at  the  zenith,  but 

an  air  mass  which  is  p  /p  times  the  vertical  one,  where  p 

*0  ^  o 

is  the  pressure  at  sea  level  and  p  is  the  pressure  at  place 
of  observation.  Therefore  the  values  which  are  obtained 
from  the  tables  do  not  refer  to  the  vertical  air  mass  of 
unity.  In  order  to  compare  data  at  stations  with  different 
elevations,  the  value  of  T  must  be  related  to  the  air  mass 
vertically  above  a  station  at  mean  sea  level.  For  stations 
which  have  a  mean  pressure  that  is  more  than  50  millibars 
lower  than  the  mean  pressure  at  sea  level,  the  value  of  T 
should  be  reduced  at  sea  level  with  the  aid  of  the  following 
relationship  given  in  [10]: 


Tp  =  1  +  (T  -  1)  •  (P(mh)/P(m)) 


(4.7) 


* 


97 


where  is  the  turbidity  factor  extrapolated  to  sea 

level , 

P (m^)  is  the  tabular  value  of  the  function  P (m)  at 
the  relative  air  mass  m^,  and 

P (m)  is  the  tabular  value  of  the  function  P (m)  at 
the  optical  air  mass  m. 

4.1.2  Angstrom's  turbidity  coefficient,  3:  a  more 
appropriate  way  of  expressing  the  scattering  of  dust  in  the 
atmosphere  independent  of  the  molecular  scattering  for  the 
non  selective  part  of  the  solar  spectrum  has  been  introduced 
by  Angstrom,  1929  [2].  In  the  ultraviolet  and  visible  portion 
of  the  solar  spectrum,  where  there  is  almost  no  selective 
absorption  apart  from  the  ozone  absorption  at  high  layers, 
attenuation  of  solar  radiation  is  caused  only  by  dry  air  and 
dust  in  the  atmosphere.  For  this  part  of  the  spectrum  the 
combined  effect  of  the  scattering  by  dry  air  molecules  and 
dust  may  be  expressed  by  the  following  equation: 

1^  =  Iox  •  exp-(aR(A)  +  aD  ( A )  )  *m  (4.8) 

where  I,,  I_, ,  a  (A)  and  m  are  as  defined  before,  and 
A  U  A  R 

aD(A)  is  the  extinction  coefficient  of  the  atmospheric 
dust. 

The  separate  effects  of  absorption  and  scattering  are 

O 

difficult  to  separate  for  dust,  but  as  a  compromise.  Angstrom 
[3]  found  that  the  extinction  coefficient  of  the  dust  in  the 
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atmosphere  may  be  approximated  in  a  form  similar  to  that  of 
the  Rayleigh  scattering  coefficient.  In  1929,  1930  he  proposed 
the  well  known  empirical  relation  for  the  wavelength  dependence 
of  extinction  of  solar  radiation  by  dust  given  by: 

aD(X)  =  6  •  X_Y  (4.9) 

where  3  is  the  turbidity  coefficient, 
y  is  a  wavelength  exponent,  and 
X  is  the  wavelength. 

The  magnitude  of  the  wavelength  exponent  y  depends  on  the 
average  size  of  the  dust  particles  in  the  atmosphere.  For 
small  particles  like  those  demanded  by  the  theory  of  Rayleigh, 
y  may  be  as  large  as  4.0.  For  larger  particles  it  may  decrease 
to  a  value  between  zero  and  1.0.  In  practice  y  remains 
between  the  limits  0.5  and  3.  Angstrom  [4]  found  that  y  = 

1.3  ±  0.2  to  be  a  reasonable  average  value. 

The  turbidity  coefficient  3  depends  upon  the  amount  of 
dust  present  in  the  atmosphere  and  for  a  constant  value  of  y, 
it  is  proportional  to  the  amount  of  dust  in  the  air.  It  is 
clear  from  equation  4.9  that  the  quantity  3  coincides  with 
the  dust  extinction  coefficient  aQ  for  wavelength  X  =  ly. 

The  method  most  widely  used  for  determination  of  the 

O 

turbidity  coefficient  is  Angstrom's  original  method  recommen¬ 
ded  in  the  IGY  Instruction  Manual  (1958) .  It  is  based  on  the 
fact  that  the  extinction  by  dust  in  the  atmosphere  is  assumed 


- 

. 
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to  take  place  according  to  the  criteria  for  equations  4.8  and 
4.9,  and  the  measurement  of  the  radiation  behind  the  Schott 
filter  of  RG2  glass  and  the  total  radiation  without  filter. 

The  RG2  filter  transmits  more  than  90  per  cent  of  all 
radiation  of  wavelengths  above  0.63  micron.  Practically  the 
total  water  vapor  absorption  in  the  solar  spectrum  falls  on 
the  wavelengths  which  are  transmitted.  If  the  total  radia¬ 
tion  and  the  radiation  behind  the  RG2  filter  (X  >  .  63y)  are 
measured,  then  the  short  wave  radiation  (.2  9  <  A  <  .  6 3 pi )  can 

be  obtained,  as  discussed  in  chapter  3,  as  the  difference 
between  the  total  radiation  and  the  filtered  radiation.  The 
effect  of  water  vapor  is  eliminated. 


Integrating  equation  4.8  for  wavelengths  between  zero 
and  0.6  3y  and  applying  the  reduction  factor  R  for  mean  sun- 
earth  distance,  the  short  wave  radiation,  Ik,  is  expressed 
by: 


0.630 


Ik  = 


VX) 


exp  (-m  aR(  A) 


-  m 


aDU)  ) 


dA 


(4.10) 


0 


The  integral  value  of  4.10  is  represented  in  tables  such  as: 

-2  -1 

Tables  for  Radiation  intensity  of  the  sun  (m  cal  .cm  .  m  ) 
in  relation  to  air  mass  m  and  the  turbidity  coefficient  3 
for  wavelength  less  than  0.630y  computed  by  Feussner  (1932) 
or  Hoelper  (1939).  The  data  in  Table  7,  given  in  the  appen¬ 
dix  to  the  IGY  Manual  [10] ,  were  calculated  from  these  sources 
with  due  regard  to  the  international  Pyrheliometric  Scale 
1956,  and  the  extraterrestrial  solar  spectrum  by  Nicolet 
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(1951) .  Also  there  are  diagrams  which  represent  the  value  of 
the  integral  4.10  for  different  air  masses  and  different 

O 

values  of  3  like  those  given  by  Angstrom  (1929)  [2].  If  the 

intensity  of  the  short-wave  radiation  at  different  values  of 
air  mass  m  corresponding  to  different  times  at  any  place  are 

O 

calculated,  then  Angstrom  turbidity  coefficient  can  be  inter¬ 
polated  from  the  standard  table  or  diagrams.  Usually  tables 
and  diagrams  for  the  turbidity  coefficient  3  are  built  up  to 
refer  to  a  place  of  observation  with  the  standard  pressure, 
at  mean  sea  level,  of  1000  millibars.  For  elevated  stations 
with  atmospheric  pressure  lower  than  950  millibars  the  value 
of  the  turbidity  coefficient  3  should  be  reduced  to  standard 
pressure  using  the  relation: 

Bp  =  3  •  p/pQ  (4.11) 

In  integral  4.10,  a  somewhat  larger  value  for  a  should 
be  used  to  account  for  ozone  absorption  in  the  ultraviolet 
and  visible  regions  of  the  solar  spectrum. 

4.1.3  SCHUEPP'S  TURBIDITY  COEFFICIENT,  B:  Angstrom's 
turbidity  coefficient  3  refers  to  a  wavelength  of  l.Oym, 
which  is  outside  the  visible  spectrum.  Schiiepp  (1949)  has 
substituted  a  turbidity  coefficient  B  which  refers  to  a 

O 

wavelength  of  0.5ym.  Angstrom  [7]  found  by  observation  that 
in  the  center  part  of  the  visible  spectrum 

B  =  0.434  3  •  2T 
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The  coefficients  B  and  8  are  not  related  by  a  constant  factor, 
but  depend  specifically  on  the  value  of  y.  However,  calcu¬ 
lations  for  B  are  often  made  by  letting  y  =  1.3  as  suggested 

O 

by  Angstrom.  Then 
B  =  1.069  6 

A  direct  measurement  of  B  can  be  obtained  by  using  the  Volz 
type  G  Sunphotometer . 


4.2  DETERMINATION  OF  LINKE'S  TURBIDITY  FACTOR  FOR  EDMONTON 

Table  6  of  the  IGY  Instruction  Manual  1958  [10]  was 
divided  into  several  intervals  with  each  interval  covering 
one  unit  of  air  mass.  A  polynomial  was  fitted  for  each 
interval  between  the  values  of  the  air  mass  m  and  the  factor 
P (m)  using  a  least  square  technique.  Table  4.1  shows  the  ten 
polynomials  used  for  calculating  the  factor  P (m)  for  any  given 
value  of  air  mass  m  in  this  study. 

The  turbidity  factor  T  for  Edmonton  was  determined  from 

-2  -1 

equation  4.6  using  the  solar  constant  Iq  =  1.98  cal  .  cm  .  min  , 

the  values  for  the  reduction  factor  R  from  reference  [25] , 

the  measured  total  direct  normal  solar  radiation  at  Edmonton 
-2  -1 

in  cal  .  cm  .  min  and  the  calculated  value  of  P (m) .  Then 
the  extrapolated  turbidity  factor  T^  was  determined  according 

to  equation  4.7.  For  each  day,  T  was  calculated  at  hourly 

P 

intervals  starting  from  hour  angle  zero,  at  solar  noon,  for 


TABLE  4.1 


Polynomials  for  calculating  the  factor 
P (m)  as  a  function  of  the  air  mass  m 


P  (m) 

2 

=  112.39  -  213.96  .  m  +  180.99  .m 
-  56.08  •  m3 

0  <  m  <  1 

P  (m) 

2 

=  55.81  -  49.78  •  m  +  20.06  .m 
-  2. 95  •  m3 

1  <  m  <  2 

P  (m) 

2 

=  25.97  -  5 . 6  5  •  m  -  1.45  •  m 

3  4  5 

+  0.024  •  m  +  0.46  •  m  —  0.131  •  m 

c 

+  0.01  •  m 

2  <  m  ^  3 

P  Cm) 

=  8.81  +  5.06  •  m  -  2.48  •  m2  +  0.29  •  m3 

3  <  m  ^  4 

P  (m) 

2  3 

=  8. 06  +  2. 52  .  m  -  1. 04  .  m  +0.09.  m 

4  <  m  <  5 

P  (m) 

2 

=  13.64  -  2.02.  m  +  0.12. m 

5  <  m  <  6 

P  (m) 

=  9.11  -  0.59  .  m  +  0.004  .  m2 

6  <  m  <  7 

P  Cm) 

=  8.031  -  0.408  •  m 

7  <  m  <  8 

P  Cm) 

=  7.309  -  0.318  •  m 

8  <  m  <  9 

P  Cm) 

=  6.752  +  0.256  •  m 

9  <  m  <  10 

. 


♦ 
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both  forenoon  and  afternoon  periods.  Figure  4.1  shows  the 
daily  average  value  of  for  Edmonton  and  the  mean  monthly 
values.  The  latter  are  listed  in  table  4.2. 


TABLE  4 . 2  Mean  monthly  value  of  Linke's 
turbidity  factor  T^  for  Edmonton 


Month 

T 

P 

Month 

Tp 

January 

2.01 

July 

3.26 

February 

2 . 08 

August 

3.22 

March 

2.49 

September 

2.75 

April 

2.92 

October 

2.29 

May 

3.26 

November 

2.09 

June 

3.37 

December 

1.98 

T  is  lowest  in  December  and  January  and  increases  to 
P 

a  flat  peak  in  the  months  May  through  September.  The  in¬ 
crease  between  winter  and  summer  values  of  T^  is  partially 
due  to  the  increase  of  absorption  in  the  near-infrared  bands. 


The  mean  yearly  value  of  T  ,  averaging  all  data  is  2.64.  The 
merit  of  T  is  that  it  can  be  found  from  measurement  of  the 


P 

intensity  of  total  solar  radiation.  However,  T^  is  affected 
by  absorption  of  solar  radiation  in  the  near-infrared  water 


. 
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Figure  4.1  Annual  variation  of  Linke's  turbidity  factor  for  Edmonton 
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vapor  bands.  Therefore,  T  can  not  essentially  be  considered 
as  a  strict  measure  of  the  atmospheric  turbidity  alone  [24] . 

4 . 3  DETERMINATION  OF  ANGSTROM'S  TURBIDITY  COEFFICIENT  FOR 

EDMONTON 

O 

Angstrom  turbidity  coefficient  3  has  been  computed  for 

O 

Edmonton  according  to  Angstrom's  original  method.  First, 

the  hourly  values  of  the  radiation  measured  by  the  RG2  filter 

were  prepared  (as  discussed  in  Chapter  3) .  Then  these  values 

were  subtracted  from  the  corresponding  hourly  values  of  the 

-2  -1 

total  normal  incidence  radiation  (in  cal.  cm  .  min.  )  and 

the  differences  were  reduced  to  mean  solar  distance  so  the 

shortwave  radiation  defined  by  X  <  630my  was  obtained.  This 

radiation  can  be  regarded  as  free  from  selective  absorption 

and  was  used  for  determination  of  3  by  interpolation  from 

table  7  given  in  the  appendix  of  the  IGY  Instruction  Manual. 

This  table  is  based  on  the  use  of  y  =  1.3.  Then  the  3  values 

thus  derived  were  reduced  to  3  according  to  equation  4.11. 

P 

O 

Table  4.3  lists  the  monthly  average  value  of  Angstrom's 

turbidity  coefficient  3  for  Edmonton  and  figure  4.2  shows 

P 

the  annual  variation  of  3  .  The  annual  trend  of  3  is 

P  P 

similar  to  that  for  Tp.  The  peak  in  3p  occurs  in  May  and 
June  reaching  an  average  value  of  0.048  and  the  lowest  value 
of  0.01  occurs  in  the  winter  months.  The  mean  value  of  3 

P 

over  the  year  was  found  to  be  0.026. 
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Figure  4.2  Annual  variation  of  Angstrom  turbidity  factor, ,  for 
Edmonton. 
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TABLE  4 . 3  Mean  monthly  value  of  Angstrom' s 


turbidity  factor  for  Edmonton 


Month 

bp 

Month 

sp 

January 

0.012 

July 

0.041 

February 

0.019 

August 

0.037 

March 

0.034 

September 

0.018 

April 

0.038 

October 

0.012 

May 

0.047 

November 

0.008 

June 

0.048 

December 

0.010 

The  seasonal  variation  of  3  with  a  pronounced  maximum 

P 

in  the  early  summer  months  is  a  well  known  phenomena  [24] . 
The  obvious  cause  of  this  must  be  that  on  account  of  the 
rapid  heating  of  the  ground  in  spring  and  early  summer  and 
that  consequently  dust  particles  are  then  carried  most 
effectively  from  the  ground  to  the  higher  layers  of  the 
atmosphere,  while  during  winter  the  surface  is  covered  by 
snow  much  of  the  time,  suppressing  the  dispersion  of  surface 
dust. 


4.4  ABSORPTION  OF  SOLAR  RADIATION  BY  ATMOSPHERIC  WATER  VAPOR 

In  the  near  infrared  regions  of  the  solar  spectrum  (A> 
630  my) ,  absorption  by  water  vapor  and  absorption  by  carbon 
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dioxide  dominate  the  attenuation  of  radiation.  This  attenua¬ 
tion  occurs  by  a  series  of  absorption  bands  each  comprised 
of  hundreds  of  narrow  individual  absorption  lines,  which 
merge  to  form  the  bands.  The  wings  of  these  numerous  lines 
contribute  a  definite  amount  to  the  continuum  extinction. 

So  the  infrared  extinction  is  comprised  of  two  parts:  a 
continuum  extinction  at  all  frequencies,  caused  by  scattering 
and  the  wings  of  absorption  lines,  and  a  selective  absorp¬ 
tion  within  each  individual  line  contributing  to  the  total 
absorption  band.  The  absorption  is  mainly  caused  by  the 
water  vapor  and  is  simply  denoted  as  water  vapor  absorption. 

Different  formulas  for  the  calculation  of  the  water 
vapor  absorption  in  the  infrared  regions  of  the  solar  spec¬ 
trum  have  been  derived  by  scientists.  As  a  compromise  the 
water  vapor  absorption  is  expressed  as  a  function  of  the 
total  amount  of  water  vapor,  in  the  zenith  direction,  betwe¬ 
en  the  surface  of  the  earth  and  space,  known  as  the  precipi- 
table  water.  The  amount  of  precipitable  water,  in  centi¬ 
meters,  is  the  thickness  of  liquid  water  which  would  be 
formed  if  all  the  water  in  the  zenith  direction,  vertical 
column,  were  condensed  at  the  surface. 

O 

In  1961  Angstrom  [5]  has  derived  the  following  formula 
for  the  infrared  absorption  (the  region  where  the  main  water 
vapor  bands  are  found  is  0.75  to  2.1  y)  by  the  water  vapor 
in  the  case  of  a  perfectly  clear  atmosphere  (turbidity 
coefficient  3=0)  at  standard  sea  level  pressure 


:'*rr 

* 


WA  =  0.31  -  0.21  exp (-0.23  w*m^) 
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(4.12) 

where  WA  is  the  infrared  absorption  by  the  water  vapor  (in 

,  -2  .  -1. 
cal .  cm.  mm .  )  , 

w  is  the  precipitable  water  content  of  the  vertical 
column  (in  cm)  and 
m^  is  the  relative  air  mass. 


In  1964  Angstrom[6]  expanded  his  formula  to  take  regard  also 
for  the  influence  of  scattering  within  the  wavelength  region 
where  the  main  water  vapor  bands  are  found,  viz.  at  0.75  to 
2.1  y.  He  also  added  two  corrections  for  the  pressure  depen¬ 
dence  of  the  absorption  and  for  absorption  above  2.1y.  He 
proposed  that  for  m^.w  between  0.5  and  10  and  6  <  1.0,  the  infra¬ 
red  absorption  by  water  vapor  may  be  calculated  by  the  equation: 

0.30  m,  •  6 

WA  =  0.163  (w-mh)  .  (0.81)  (4.13) 

O 

where  6  is  Angstrom's  turbidity  coefficient.  This  equation  is 
applicable  to  the  average  conditions  at  stations  near  to  sea 
level  with  uncertainty  in  deriving  values  for  WA  of  at  least 
±5  percent  [6] . 

4 . 5  DETERMINATION  OF  WATER  VAPOR  ABSORPTION  AT  EDMONTON 

The  water  vapor  absorption,  WA,  at  Edmonton  was  calcula¬ 
ted  according  to  equation  4.13,  using  the  calculated  relative 


.  IS, 
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air  mass  values  for  Edmonton,  the  average  daily  precipitable 

water  values  for  forenoon  and  afternoon  obtained  from  reference 

[11]  and  the  3  values  determined  for  Edmonton  in  section  4.4. 

For  each  day  WA  was  calculated  at  hourly  intervals  for  both 

forenoon  and  afternoon  periods  then  the  daily  average  value 

was  calculated.  Table  4.4  lists  the  monthly  mean  value  of  the 

water  vapor  absorption,  WA,  at  Edmonton  and  figure  4.3  shows 

the  annual  variation  of  WA.  The  atmospheric  water  vapor 

absorption  shows  higher  values  in  the  summer  months  with  the 

-2  —1  —2  -1 

maximum  of  0.245  cal. cm.  min  (54.19  Btu.ft.  hr  )  in 

August.  The  lowest  values  occur  in  the  months  October  to 

-2  .  -1 

March  with  an  approximate  average  value  of  0.18  cal. cm.  min 
(39.82  Btu.  ft-2  .hr-JL)  .  These  correspond  to  an  amount  of 
precipitable  water  of  1.84  cm  (0.724  inches)  and  0.48  cm 
(0.19  inches),  respectively. 


V 


TABLE  4 . 4  Mean  monthly  value  of  water  vapor 


absorption  at  Edmonton 


VIA 

Month 

n  -2.-1 
cal .  cm  .  mm 

-2  -1 
Btu.ft.  hr 

January 

0.178 

39.37 

February 

0.180 

39.82 

March 

0.181 

40.03 

April 

0.182 

40.26 

May 

0.185 

40.92 

June 

0.210 

46.45 

July 

0.215 

47.56 

August 

0.243 

53.75 

September 

0.218 

48.22 

October 

0.184 

40.70 

November 

0.180 

49.82 

December 

0.176 

38.93 

V 

WA  (Btu  ft 
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Figure  4.3  Annual  variation  of  water  vapor  absorption,  WA  (Btu  ft  hr  ), 
at  Edmonton. 


CHAPTER  5 


DISCUSSION  AND  CONCLUSIONS 
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In  this  study,  solar  radiation  measurements  at  Edmonton, 
Alberta,  Canada  were  analyzed  for  the  period  from  May  1972  to 
September  1975.  The  data  consisted  of  records  of  clear  sky 
direct  normal  incidence  radiation  within  the  whole  spectrum 
and  its  spectral  components  passing  through  the  RG2  and  RG8 
filters.  The  spectrum  was  divided  into  five  bands  defined 
by  the  wavelength  intervals:  A  <  626my  (shortwave),  A  >  626my 
(red  and  infrared) ,  A  <  690my  (visible) ,  A  >  690my  (infrared) 
and  626  <  A  <  690my  (red) .  In  order  to  study  the  effect  of 
the  atmospheric  scattering  and  water  vapor  absorption  on  the 
solar  energy  more  attention  was  given  here  to  the  shortwave 
and  infrared  bands.  The  intensity  of  the  shortwave  radiation 
received  at  the  earth's  surface  is  mainly  dependent  on  atmos¬ 
pheric  scattering  while  the  intensity  of  the  infrared  radiation 
is  mainly  dependent  on  the  absorption  by  the  atmospheric  water 
vapor . 

5.1  APPARENT  SOLAR  CONSTANT  AND  APPARENT  EXTINCTION  COEFFICIENT 

* 

The  apparent  solar  constant  I  and  the  apparent  extinc¬ 
tion  coefficient  A  have  been  calculated  from  the  total  direct 
normal  incidence  radiation  IDN  and  the  optical  air  mass  m  using 
the  equation: 


* 

IDN  =  I  exp (-Am) 

* 

The  annual  variation  and  the  monthly  average  values  for  I  and 


' 


■ 


115 


A  are  reported  in  section  2.4.  The  mean  annual  values  obtained 

-2  -1  * 

are  346  Btu.ft  .hr  .  for  I  and  0.122  for  A.  The  results  are 

compared  with  the  values  obtained  by  Sadler  [28]  for  the 

same  location  for  the  period  from  1965  to  1968.  The  present 

* 

results  show  about  5  percent  lower  values  for  I  for  April  to 

August,  inclusive.  A  better  agreement  appears  for  the  rest 

of  the  year.  The  present  values  for  A  are  slightly  lower  for 

January  to  September  and  slightly  higher  for  October  to 

December.  The  simplicity  of  this  method  used  to  represent  the 

solar  insolation  is  desirable  in  practice,  but  additional 

investigation  is  required  to  determine  the  effect  of  the 

* 

climatic  conditions  on  the  values  of  I  and  A.  Since  the 

recording  station  is  located  near  the  center  of  the  city, 

consideration  should  be  given  to  the  dust  content  in  the 

atmosphere  in  order  to  arrive  at  more  representative  values  of 
* 

I  and  A. 


5.2  DIURNAL  AND  ANNUAL  VARIATIONS  OF  THE  DIRECT  RADIATION 
AND  ITS  SPECTRAL  COMPOSITION 

The  annual  variation  of  the  total  direct  normal  incidence 
radiation  and  its  composition  within  the  five  spectral  bands 
at  hour  angles  0  to  6  is  shown  graphically  in  sections  2.6  and 
3.5.  The  monthly  average  values  are  also  reported.  These 
values  occur  at  about  the  third  week  of  the  month  but  measure¬ 
ments  on  adjacent  days  may  show  up  to  10  percent  variation 
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from  the  mean.  The  maximum  intensity  of  the  total  direct 

-2  -1 

radiation  of  about  305  Btu.ft  .hr  occurs  in  March  and  April 

at  solar  noon.  The  maximum  intensity  of  the  infrared  radia- 

-2  -1 

tion  of  about  182  Btu.ft  .hr  occurs  also  in  March  and  April 
when  the  precipitable  water  in  the  atmosphere  is  low.  Lower 
energy  values  are  found  in  the  summer  months  as  the  precipi¬ 
table  water  level  increases  in  the  atmosphere.  The  maximum 

-2  -1 

intensity  of  the  shortwave  radiation  of  about  125  Btu.ft  .hr 
occurs  in  the  summer  months  at  solar  noon  when  the  path  length 
through  the  atmosphere  is  relatively  short.  The  available 
energy  decreases  as  the  hour  angle  deviates  from  solar  noon. 

The  decrease  is  more  rapid  after  hour  angle  2  due  to  rapid 
increase  in  the  path  length  of  the  radiation  through  the 
atmosphere . 

The  infrared  fraction  of  the  direct  radiation  has  minimum 
value  during  the  summer  months  (for  all  hour  angles) ,  although 
the  path  length  through  atmosphere  is  relatively  short.  This 
is  due  to  the  larger  quantities  of  precipitable  water  during 
this  season.  The  red  fraction  of  the  direct  radiation  has  a 
mean  value  of  6.5  percent  with  slight  annual  variation  at  all 
hour  angles.  The  shortwave  fraction  of  the  direct  radiation 
has  a  maximum  during  the  summer  months  and  decreases  gradually 
during  the  rest  of  the  year,  a  characteristic  which  occurs  at 
all  hour  angles.  This  fraction  decreases  as  the  hour  angle 
increases  from  solar  noon  because  of  the  increase  in  the  path 
length  through  the  atmosphere.  As  a  result,  the  infrared 
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fraction  increases  approximately  in  proportion  to  the  decrease 
in  the  shortwave  fraction  with  increasing  hour  angle. 


5 . 3  VARIATION  OF  THE  INFRARED  AND  SHORTWAVE  FRACTIONS  WITH 

SOLAR  ALTITUDE  AND  OPTICAL  AIR  MASS 

The  fraction  of  the  direct  radiation  within  the  infrared 
band  increases  from  about  54  to  62  percent  while  the  fraction 
in  the  shortwave  band  decreases  from  about  40  to  35  percent 
as  the  solar  altitude  decreases  from  60  to  20  degrees.  For 
further  decrease  in  altitude,  these  values  change  faster 
reaching  74  percent  within  the  infrared  band  and  20  percent 
within  the  shortwave  band  at  5  degrees. 

By  representing  the  infrared  and  the  shortwave  fractions 
as  functions  of  the  optical  air  mass,  linear  relationships  seemed 
reasonable  and  have  been  computed  by  the  least  squares  method. 

For  the  same  optical  air  mass,  the  infrared  fraction  is  less  in 
the  summer  than  in  the  winter  by  about  2.5  percent.  On  the 
other  hand,  the  shortwave  fraction  is  about  1.5  percent  higher 
in  the  summer. 


5 . 4  ATMOSPHERIC  TRANSMISSION  FACTORS  AND  EXTINCTION 
COEFFICIENTS 


For  the  same  solar  altitude  and  approximately  the  same 
optical  air  mass,  the  intensities  of  direct  radiation  within 


118 


the  whole  spectrum  and  the  five  bands  reduced  to  mean  earth- 
sun  distance  are  higher  in  winter  than  in  summer.  This  is 
consistent  with  higher  atmospheric  transmission  factors  in 
winter  when  the  precipitable  water  and  the  dust  content  in 
the  atmosphere  are  lower  than  in  summer.  The  ratio  of  the 
extinction  coefficient  in  summer  to  that  in  winter  is  1.29 
for  the  whole  spectrum,  1.47  for  the  infrared  band  and  1.16 
for  the  short  band.  In  summer,  the  ratio  of  the  extinction 
coefficients  of  the  infrared  band  and  the  shortwave  band  to 
the  extinction  coefficient  for  the  whole  spectrum  are  0.79  and 
1.15,  respectively.  The  corresponding  values  in  winter  are  0.66 
and  1.44,  respectively. 

5.5  ATMOSPHERIC  SCATTERING  AND  ABSORPTION 

To  render  the  results  from  Edmonton  directly  comparable 
with  those  from  other  places,  it  was  regarded  appropriate  to 
determine  some  parameters  that  facilitate  the  treatment  and 
give  the  results  a  general  character.  Three  such  parameters 

O 

are:  Linke's  turbidity  factor.  Angstrom's  turbidity  coefficient 

and  the  water  vapor  absorption. 

Linke's  turbidity  factor  T^  has  shown  a  marked  annual 

variation  with  high  values  in  summer  and  low  values  in  winter. 

The  mean  value  over  the  year  is  2.64.  The  increase  between 

winter  and  summer  values  of  T  is  due  to  the  increase  of  dust 

P 

and  other  aerosol  particles  in  the  atmosphere.  It  is  also 
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partially  due  to  the  increase  of  atmospheric  water  vapor  in 
the  atmosphere  in  the  summer. 

O 

The  annual  trend  of  Angstrom's  turbidity  coefficient  3^ 

is  very  similar  to  that  for  T  ,  and  for  the  same  reasons.  The 

P 

yearly  mean  value  for  3  is  0.026. 

P 

The  atmospheric  water  vapor  absorption  is  found  to  be 
-2  -1  -2  -1 

47.0  Btu.ft  .hr  (0.213  cal  cm  min  )  with  an  annual  varia¬ 
tion  indicating  a  maximum  in  July  and  August  and  a  minimum  in 
the  winter  months.  This  corresponds  to  an  amount  of  precipitable 
water  of  0.465  inches  (1.16  cm)  in  the  vertical  column  of  the 
atmosphere . 


* 
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